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Outline  of  the  course  w«  t  as  /o  !L,  s  ; 


D.  Stanat  presented  an  introduction  to  algorithms  and  data  structures  of  general  interest 

^  () 

T.  Ottmam  introduced  the  fleld  of/icomputatlonal  geometryjwith  a  talk  on  geometric  algorithms  in  the 
plane.  He  demonstrated  geometrical  dlvlde-and<onquer  ^  the  prdblem  of  reporting  intersecting 
linesegments.  Then  he  discussed  sweep  tine  algorlthms^srme  same  problems.  Introducing  the  use  of 
skeleton  structures  and  exploring  minlmixation  op;iiemdry  space  requirements.  The  usage  of  plane  sweep 
to  report  intersecting  Iso-orientedrectgDglea-letrto  the  introduction  of  segment  trees,  interval  trees  and 
priority  search  trBa^lflsao  dwaH^^tlon  of  a  graph  was  demonstrated  -on  Uie  problem  of  polygon 
intersection.  FlnMIyia  hidden  line  detection  dgorlthm  was  presented. 

J.  Hopcroft  discussed^eometrical  problems  related  to  robotics  After  giving  an  overview  of  the  field,  his 
talks  centered  on  the  two  problems  of  automated  generationo^endlng  surfaces  and  on  motion  planning. 

P.  Widmaytr  discussed^ieurisUcs  for  flrtding  approximations  for  Steiner  minimum  treesj  The  exact 
solution  of  this  pro^lm  Is  kQown  to  be  np-complete.  J 

K,  Hlnrichs  descrlbed|the  grid  file  as  a  data  structure  suited  for  geometrical  computation' and  presented 

experience  with  the  Imidementation.  ^ 

cr^p- - — - 

A.  UeUr  ex^alned/tdifferent  schemes  for  representing  three-dimensional  object^  Furthermore  he 
discussed  how  the  Relational  Data  Base  model  could  be  used  to  store  geometric  objects.  He  pointed  out 
some  of  the  problems  of  this  approach  and  indicated  possible  extensions  of  ther^attoif^  model  to  make  it 
suitable  for  computational  geometry.  _ _ _ _ — ^ 


F.  Luccio  discussed  ^visibility  problems  that  occur  in  VLSI  deslgi^' He  showed  the  equivalence  of  a 
visibility  problem  with  a  planar  graph  and  hence  the  appllcabijl^  of  graph  theory.  He  also  briefly 
explain^  the  view  of  a  program  as  a  decision  tree.  ^ 

7) 

F.  Prepamta  In  his  overview  of  geometric  algorithms  introduced^ugorithms  for  point  location,  convex  hull 
and  the  maxima  of  a  set  of  vectors  in  two  and  three  dimensions.  .[Die  time  complexity  of  these  algorithms 
was  discussed.  Algorithms  for  constructing  Voronoi  diagranmik  were  introduced  to  solve  proximity 
problems  and  the  use  of  plane  sweep  algorithms  for  solving  planar  intosectlon  problems  was  shown. 

^ - - 

Additional  contributions  c^me  flrom  J.  Sack  and  Th.  Strothotte  who  presented  a  recent  result  ommerging 
heaps  and  also  pointed  ouVsome  unsolved  problems.  ^ 

Further  information  on  the  talks  can  be  found  in  the  enclosed  papers  handed  out  by  J.  Hopcroft,  K. 
Hlnrichs,  A.  Meier  and  F.  Luccio. 
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Prof.  J.  Nleveigelt  f 


Copy 


.  j  .1-  : 


^>1?-  ^  .  ’ 


•2-  ^/^(rncz-j, 


t-33100  UDINE  (Italy).  Palazzo  dal  Torso.  Pl.izza  Garibaldi,  18 
Tel.  (04321  29  49  69  -2  25  23 


Course  on 


ALGORITHMS  AND  DATA  STRUCTURES  FOR  GEOMETRIC  COMPUTATION 

Udine.  July  8-12,  1985 

TIME  TABLE 


9.00  -  10.  30  Registration 

10.  30  -  12.  30  D.  Stanat:  Introduction  to  Algorithms  and 

Data  Structures,  I  ; 

14.  30  -  16.  30  D.  Stanat:  Introduction  to  Algorithms  and 

Data  Structures,  II  ; 


Tuesday:  9.  00  -  1  1 .  00  T.  Ottmann:  Geometric  Algorithms  in  the 

Plane,  1  ; 

1  1.  30  -  12.  30  F.  Luccio:  Visibility  Problems,  I; 

14.  30  -  16.  30  J.  Hopcroft;  Robotics  Algorithms,  I; 

Wednesday:  9.00  -  11.00  T.  Ottmann:  Geometric  Algorithms  in  the 

Plane.  II; 


Thu  rsda^ 


F  ridaT 


1 1.  30  -  12.  30 


K.  Hinrichs  :  Data  Structures  for 
Geometric  Computation; 


14.  30  -  16.  30  J.  Hopcroft:  Robotics  Algorithms,  II; 

^  9.  00  -  11. 00  A.  Meier;  Data  Bases  for  Geometric  Objects; 

1  1 .  30  -  12.  30  F.  Luccio:  Visibility  Problems,  II; 

14.  30  -  16.  30  F.  Preparata:  Overview  on  Geometric 

Algorithms,  I; 

9.  00  -  11.  00  F.  Preparata:  Overview  on  Geometric 

Algorithms,  II; 

11.  30  -  12.  30  Discussions. 

Possible  contributed  talks  will  be  delivered  at  the  end  of  the 
afternoon  sessions. 
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tract  for  talks  at  CISM 

Introduction  to  Algorithms  and  Data  Structures 
Donald  F.  Stanat 

The  University  of  North  Carolina  at  Chapel  Hill 
Chapel  Hill.  North  Carolina.  USA 

se  talks  will  provide  a  fast-paced  introduction  to  the  field  of 
orithms  and  data  structures.  We  will  begin  with  a  description  of 
basic  abstract  data  structures:  arrays,  lists,  trees  and  graphs, 
ether  with  implementations  of  the  structures.  We  will  include 
sures  of  cost,  both  in  time  and  space,  of  various  operations  on 
se  structures  for  each  of  the  implementations.  Finally,  we  will 
cuss  some  selected  advanced  data  structures,  including  hash  tables, 
ps  and  a  number  of  different  kinds  of  balanced  trees. 

^  we  will  survey  some  general  algorithm  types,  including  greedy 
orithms,  divide  and  conquer,  and  exhaustive  search,  including  some 
the  ways  of  making  exhaustive  search  feasible:  backtracking,  branch 
bound,  and  dynamic  programming. 

ally  we  will  describe  the  cost  of  solving  really  hard  problems  and 
reduce  the  notion  of  NP  (Nondeterministic  Polynomial)  difficulty, 
time  permits,  we'll  also  talk  about  approximation  algorithms  for 
blems  whose  exact  solutions  have  NP  cost;  these  approxim.ations 
orithms  provide  a  means  of  getting  a  non-optimal  solution  at 
erate  cost  for  a  problem  whose  optimal  solution  is  too  expensive. 


Storage  and  access  structures  for  geometric  data  bases 

J.  Nlerersclt  tod  L  Hlnrtchj 
Inronnatli.  ETIi,  CH  S092  ZUrich 


Abstract 

Ccometrlc  compulation  and  data  bases,  two  hitherto  unrelated  computing  technologies,  have  begun  to 
Ir.nuencc  each  ether  In  response  to  the  growing  use  of  graphics  and  computer-aided  design.  CAD  Imposes 
a  rew  challenge  to  data  base  implementers.  A  data  base  Sj-stem  for  CAD  must  manage  "In  designer  real 
time"  large  collections  not  of  points,  but  of  spatial  objects.  In  such  a  way  that  proximity  queries  (such  as 
Intersection,  contact,  minimal  tolerances)  are  arswered  effldentiy.  There  are  many  techniques  for  reducing 
the  problem  of  storing  spatial  objects  to  storing  (seta  oO  points.  Common  to  ail  of  them  Is  the  problem 
that  simple  queries  on  objects  turn  Into  complex  queries  on  points  -  much  more  complex  than  orthogonal 
range  queries. 

We  describe  In  detail  a  technique  which  Is  particularly  suitable  for  storing  geometric  objects  built  up  from 
simple  primitives,  as  they  commorJy  occur  In  CAD.  Proximity  queries  are  handled  efficiently  as  part  of 
the  accessing  mechanism  to  disk.  This  technique  is  based  on  a  transformation  of  spatial  objects  Into  points 
in  hjg.ber  -ditnenaicnaJ  parameter  spaces,  and  on  the  data  structure  grid  file  that  answers  region  queries  of 
complex  shape  efficiently.  The  grid  file  Is  designed  to  store  hlgly  dynamic  sets  of  multi -dimensional  data 
in  such  a  way  that  common  queries  are  answered  using  few  disk  accesses:  a  point  query  requires  two  disk 
accesses,  a  region  query  requires  at  most  two  disk  accesses  per  data  bucket  retrieved.  We  describe  a 
soifware  pacl.age  that  implements  the  grid  file  and  some  of  Its  applications. 

Contents 

1  Geometric  computation  and  data  bases 

1.1  Three  generations  of  computing  applications 

1.2  Ccm.putalicnal  geometry  -  theory  and  practice 

1.3  Tl'.e  conventioruii  data  base  approach  to  "non-standard"  data 

1.4  Geometric  modciirg  separated  from  storage  considerations 

2  An  ippro<ich  to  combined  geometric  modeling  and  stxarirg: 

Approximation,  trars formation  to  parameter  space,  grid  file 

2.1  Transformation  to  parameter  space 

2.2  Irtersection  qu'.-ries  lead  to  cone-shaped  search  regions 

2.3  Iv.  jJ uat; rg  region  queries  with  a  grid  file 

2.4  TTie  grid  file  software  package 
2J  Cc-e  studies  of  applications 

1  Geometric  computation  and  data  bases 

Gfometric  lompjtatlon  and  data  ba.ses,  two  hitherto  unrelated  computing  technologies,  have  begun  to 
IruPuence  each  other  In  response  to  the  growing  use  of  graphics  and  computer-aided  design  (CAD).  We 
recall  their  origins,  goals,  typical  technlque.s,  and  explain  the  difficulties  each  of  them  has  in  handling  the 
requirerr.ents  of  the  other. 

1.1  Three  jjcncratlotis  of  computinc  applications 

Th.e  types  of  cemputer  applications  dominant  at  different  times  may  be  classified  Into  three  gencTatlOTiS 
a-'ccrdlrg  to  their  influence  on  the  development  of  computing. 


ProrMdings  of  the  Internfltifmei  Oinb-renee 
on  Foundatiuns  of  I)aUi  (Irgani/atifin 


May  tll-iit.  I'th.).  Kvotn,  . Japan 


L  USING  A  RELATIONAL  DBMS  FOR  SOLID  MODELING 


rraditionally,  engineering  and  design  data  has  been  handled  by  ad-hoc  or  simple  file 
ystems  with  the  inherent  disadvantages  of  high  redundancy  and  poor  or  non-existent  data 
ndependence.  In  other  words,  the  way  the  data  is  physically  organized  within  the  file  must 
)e  known  to  the  programs  that  access  the  data:  Any  change  in  the  data  organization  requires 
■hanging  the  programs  and  vice-versa.  When  the  number  of  files  increases,  the  consistent 
reatment  of  data  becomes  a  problem  in  itself.  It  is  not  surprising  therefore  that  today’s 
levelopers  of  CAD  systems  begin  to  realize  the  importance  of  independent  data 
jrganization  [Ullraan  1982].  They  start  experiments  with  databases  although  engineering 
ipplications  exhibit  characteristics  which  impose  specific  requirements  on  existing  DBMS. 
:n  this  section,  we  discuss  how  solids  either  by  the  CSG-  or  BR-approach  may  be  stored  in  a 
relational  database  and  list  the  main  advantages  and  disadvantages. 


2.1.  Constmctiye  Solid  Geometry 

Although  none  of  the  existing  solid  representation  schemes  is  suitable  for  all  applications, 
the  CSG-scheme  provides  a  concise  way  to  store  a  volumetric  object  Halfspaces  may  be 
used  as  primitives  at  the  lowest  level.  However,  the  resulting  object  representations  are  not 
necessarily  regular  sets  [Requicha  1980]  because  of  the  unboundedness  of  primitive 
halfspaces.  Instead,  cubes,  cones,  cylinders  etc.  are  usually  used  as  primitives.  A  3D  object 
can  then  be  described  with  the  following  grammar: 


<03JECT>  ::=  <PRIMITIVE>  I 

<OBJECr>  <MOTION>  ARGUMENT  I 
<OBJECT>  <OPERATION>  <OBJECr> 
<PR!.MITIVE>  =  CUBE  I  CYUNDER  I  CONE  I  SPHERE  I  ... 


<MOTION>  : :  =  TRANSLATE  I  ROTATE  I  SCALE 

<OPERATION>  ::  =  UNION  I  INTEP.SECriON  I  DIFFERENCE 


The  semantics  of  a  CSG -representation  is  clear:  Each  subtree  represents  a  solid,  i.e.  a 
regular  set,  resulting  from  the  combinatorial  or  motion  operators  to  the  subparts.  The 
dynamic  behavior  of  the  data  structure  results  from  the  recursiveness  embedded  in  the 
grammar  rules.  In  the  relational  model,  recursiveness  has  to  be  broken  down,  and  a  solid 
may  be  described  by  several  relations:  It  has  to  be  treated  as  a  whole  at  a  high  object  level 
while  providing  its  individual  details  at  lower  part  levels. 

The  conceptual  scheme  of  the  CSG-approach  is  given  in  Fig.  2:  Each  OBJECT  in  a 
CSG -representation  consists  of  several  PARTS.  In  the  relational  model,  this  hierarchical 
structure  has  to  be  described  by  at  least  two  independent  relations  in  order  to  limit  data 
redundancy.  Furihermore,  two  generic  structures  are  imposed  by  the  CSG-scheme  [Lee  and 
Fu  1983]  which  cannot  be  defined  directly  by  relations.  First,  each  part  may  either  be  a 
TRANSFORMED-,  or  a  PRIMITIVE-,  or  a  COMBINED-PART  according  to  the  grammar 
rules,  for  instance,  a  CO.MBI NED- PART  consists  of  two  parts,  namely  the  FIRST  one  and 
the  SECOND  one  plus  the  corresponding  Boolean  OP^ATION  union,  intersection,  or 
difference.  Second,  each  primitive  part  is  either  a  CUBE  or  a  CONE  or  a  CYLINDER  etc. 


orientability.  i.e.  faces  may  intersect  only  at  common  edges  or  vertices,  and  each  edge  is 
shared  by  exactly  two  faces  etc.  Mathematically,  the  surface  of  a  solid  described  in  boundary 
representation  may  be  treated  as  a  manifold. 

After  having  discussed  the  main  schemes  for  representing  solids  we  ask  which  ones  are 
suitable  for  database  techniques; 


Fig.  1:  Using  Database  Techniques  for  Solid  Modeling. 


In  Primitive  Instancing  it  is  obvioiLs  that  each  instance  may  be  considered  as  a  record  or  tuple 
in  a  Data  Base  Management  System  (DBMS).  Since  a  shape  type  and  a  limited  set  of 
parameter  values  specify  an  object,  parametrization  does  not  involve  much  work  for 
geometric  and  topological  computation.  Therefore,  every  commercial  DBMS  might  be  good 
enough  for  describing  and  storing  a  part  family.  On  the  other  hand,  both  Spatial  Enumeration 
and  Cell  Decomposition  schemes  are  not  adequate  using  database  techniques,  especially  if  the 
objects  are  described  in  fine  resolution  or  by  a  large  number  of  cells.  The  cost  of  database 
interaction  for  object  manipulation  becomes  unreasonably  high.  Storing  objects  described  by 
Constructive  Solid  Geometry  or  Boundary  Representation  schemes  in  a  database,  however,  seems 
promising. 

This  paper  will  concentrate  on  database  aspects  for  solid  modeling.  Section  2  describes  how 
objects  given  in  CSG-  or  BR- representation  may  be  mapped  into  a  relational  database 
scheme.  In  section  3,  a  surrogate  model  is  Introduced  to  better  support  geometric  and 
topological  information.  A  proposal  for  direa  handling  of  veaors,  matrices,  and  tensors  in 
the  relational  model  is  outlined  in  section  4.  First  experiments  and  conclusions  are  given  in 
section  5. 


1.  REPRESENTATION  SCHEMES  FOR  SOLIDS 


For  Computer-Aided  Design  (CAD)  of  three-dimensional  objects  (3D  objects  or  solids),  the 
geometric  and  topological  aspects  of  part  and  assembly  specification  is  important.  In 
[Requicha  1980],  several  representation  schemes  for  solids  are  discussed  some  of  which  we 
briefly  describe. 

Primitive  Instancing  is  based  on  a  family  or  group  of  objects  where  each  member  is 
distinguishable  by  a  few  parameters.  For  instance,  the  family  of  cog  wheels  may  be 
described  by  a  type  code,  the  wheel’s  diameter  and  the  number  of  equally  spaced  cogs. 
Other  properties  of  the  objects  are  not  specified  explicitly;  they  either  are  constant 
throughout  the  family  or  they  depend  on  specified  parameters.  Primitive  instancing  lacks  the 
possibility  of  combining  representations  in  order  to  create  new  or  more  complex  schemes.  It 
is  also  difficult  or  even  impossible  to  derive  geometric  and  topological  properties  directly 
from  such  schemes.  In  practice  however,  parametrization  is  applicable  (and  still  widely  used) 
as  long  as  the  catalog  of  parameters  does  not  become  to  large. 

Spatial  Enumeration  denotes  a  scheme  where  the  embedding  space  is  divided  into  a  grid  of 
volume  elements,  and  a  solid  is  represented  by  a  list  of  occupied  grid  blocks  or  elements. 
Recently,  the  octree  encoding  (e.g.  [Meagher  1982])  as  a  hierarchical  spatial  enumeration  has 
been  discussed  as  a  representation  scheme  for  solid  modeling.  It  divides  the  space  occupied 
by  a  solid  into  eight  cubic  parts  recursively  until  a  fixed  maximal  resolution  is  reached. 
There  are  some  advantages  to  this  data  structure;  e.g.  Boolean  operations,  hidden  surface 
removal  or  interference  detection  show  linear  growth  because  all  objects  are  kept  spatially 
pre- sorted  at  all  times.  However,  when  moving  objects  are  taken  into  account,  more 
computation  is  involved. 

Cell  Decomposition  methods  are  based  on  the  results  of  triangulation  theory.  A  solid  or 
polyhedron  is  decomposed  into  disjoint  parts  of  different  dimensions.  Therefore,  operations 
and  calculations  become  easier  due  to  disjointness.  Cell  decomposition  may  be  considered  as 
a  generalized  spatial  occupancy  enumeration  where  cells  neither  have  to  lie  on  a  fixed  grid 
nor  have  a  pre-specified  size  and  shape.  In  [Bieri  and  Nef  1982],  a  recursive  sweep-plane 
algorithm  is  presented  that  enumerates  the  cells  of  all  dimensions  into  which  space  may  be 
partitioned  by  a  finite  set  of  hyperplanes.  The  described  method  is  also  suitable  to  compute 
the  Euler  characteristic,  the  volume  or  other  integral  parts  of  polyhedrons  represented  in 
Boolean  form.  Local  information  is  collected  at  every  vertex  and  summed  up  for  the  result 

Constructive  Solid  Geometry  (CSG)  denotes  a  family  of  representing  schemes  where  each 
object  is  described  as  Boolean  construction  or  combination  of  solid  components  via  the 
regularized  set  operations  such  as  union,  intersection,  and  difference.  Regularity  provides  a 
natural  formalization  of  dimension  preserving  properties  [Tilove  1980],  Le.  the  result  of  a 
Boolean  operation  of  two  solids  is  volumetric;  dangling  edges  and  faces  or  isolated  points 
are  not  allowed.  It  is  important  to  note  that  each  CSG-scheme  may  be  described  as  a  tree 
where  non- terminal  nodes  represent  operations  both  for  construction  and  transformation, 
and  terminal  nodes  denote  primitives  or  arguments  of  motion  respectively. 

Boundary  Representation  (BR)  describes  a  solid  by  its  bounding  surface  which  often  is 
subdivided  into  curvature-continous  regions  known  as  faces.  Each  face  as  a  region  of  its 
underlying  surface  is  again  bounded  by  a  perimeter  ring  of  edges  which  are,  in  turn, 
bounded  by  a  pair  of  vertices.  The  bounding  surface  has  some  unique  characteristics  such  as 
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Portability  has  been  achieved  by  deflnlr^  two  Interfaces:  One  towards  the  host  (hardware  and  operating 
system),  the  other  towards  client  programi.  The  software  can  be  transferred  to  other  systems  by  adapting  a 
module  GFHost  that  Isolates  the  machine-  and  disk-dependent  parts  of  the  program.  It  provides  procedures 
for 

-  creating  and  Initializing  disk  storage; 

-  opening  and  closing  communication  channels  between  the  disk  and  the  grid  file  module; 

-  creating,  deleting,  reading  and  writing  disk  blocks; 

-  managing  empty  disk  blocks. 

The  interface  towards  client  programs  consists  of  several  modules  that  provide  utility  and  query 
procedures; 

-  creating,  deleting,  opening  atxl  closing  a  grid  file. 

-  Inserting  and  deleting  records  in  a  grid  file. 

-  changing  non-key  Information  In  a  record. 

-  point  query:  find  all  records  with  given  key  values  xj. (If  keys  are  unique  at  moat  one  record  will 
be  found). 

-  range  query:  find  all  records  whose  key  values  xj  He  In  given  Intervals  Hi.  ui)  (1^1^  k). 

-  user  defined  region  query;  the  user  has  to  write  a  procedure  which  Is  called  by  the  grid  /lie  system  and 
determines  whether  a  grid  cell  (given  by  Intervals  [Ij.  ui)  (1  ^  I  ^  k))  Intersects  the  search  region  defined 
by  the  user. 

-  nextabove,  nextbelow:  given  key  1  with  key  value  xj,  find  the  records  with  key  values  above  or  below  xj 
and  next  to  xj;  this  gives  the  user  the  possibility  to  pnxess  the  records  sequentially  with  respect  to  one 
key. 

-  join  query:  the  join  query  Is  a  generalization  of  the  join  operator  known  from  relational  data  bases.  The 
user  has  to  write  some  procedures  which  are  called  by  the  grid  /lie  system  and  guide  the  join  query. 

-  counting;  the  above  queries  can  be  performed  by  only  counting  the  records,  but  not  transferring  them  to 
the  user  program. 

15  Case  studies  of  appUcations 

The  grid  file  software  package  has  been  used  to  store  and  process  geometric  objects  In  the  following 
applications  iHln  85al. 

Producing  layout  masks /or  Integrated  circuits  (Fig.  2.7).  Mask  generated  by  a  CAD  system  for  chip  design 
(David  .Mann  Forma!)  are  presented  as  a  set  of  aligned  rectangles.  Fabrication  requires  thal  a  mask  la 
represented  as  the  set  o/  connected  components  generated  by  rectangle  overlap,  I.  e.  a  set  of  aligned 
polygons  (all  edges  p>arallel  to  the  coordinate  axes,  Manhattan  geometry).  This  transformation  program  was 
implemented  by  pnxresslng  rectangles  In  a  4-cllmenslona]  grid  file  and  computing  the  connected 
components  by  Intersection  queries. 

Preprocessing  plotter  files.  In  a  CAD  system  for  mechanical  engineering  plotter  files  are  preprocessed  In 
order  to  reduce  ihe  total  distance  along  which  the  raised  pen  has  to  be  moved.  The  task  of  finding  an 
optimal  solution  to  this  problem  Is  equivalent  to  the  traveling  salesman  problem  and  therefore 
NT -complete.  The  ploUer  files  contain  line  segments  and  arcs  which  have  to  be  drawn.  The  end  points  of 
the  line  segrr.ents  and  the  arcs  arc  stored  in  a  2-dlmenslonaJ  grid  file.  A  reduction  of  the  total  pen  plotting 
time  Is  achieved  by  nearest  neighbor  queries  on  this  grid  file.  A  similar  method  using  quad  trees  Is 
presented  In  [And  83). 

Analyzing  phctcgrapf'jc  satelUle  data.  A  photograph  obtained  by  a  satellite  consists  of  512  •  512  plxelv 
Each  pixel  Is  assigned  four  color  values  In  the  range  from  0  to  255.  These  pixels  are  stored  In  a 
4 -dimensional  grid  file.  The  ground  imaged  by  these  pixels  Is  then  classified  Into  water,  forest,  fields, 
residential  and  metropolitan  areas  etc.,  by  range  queries  on  this  grid  file. 

Managing  simple  spatial  objects.  An  interactive  program  manages  large  sets  of  simple  spatial  objects,  e.  g. 
reclargles,  circles  and  segments.  These  objects,  each  of  which  Is  defined  by  a  fixed  number  of  parameters, 
are  stored  In  different  grid  files,  one  for  each  type  of  object.  The  p>rogram  allows  the  user  to  Insert  and 
delete  simple  spatial  objects  and  to  perform  proximity  queries  (e.  g.  Intersection,  containment)  on  the 
stored  data. 

Processing  geographic  data.  The  Swiss  Federal  Office  for  Statistics  made  available  to  us  a  file  which 
contains  raster  Information  about  Switzerland.  Each  record  In  this  file  represents  a  square  of  100  meters  by 
lOQ  meters  (1  hectare).  Switzerland  covers  about  4’000’000  hectares,  but  only  about  lOfrOGO  hectares  have 


ZJ  Eraluatins  region  qaerics  with  a  grid  fUe 

Wt  have  seen  that  proximity  queries  on  spatial  objects  lead  to  search  regions  significantly  more  complex 
than  orthogonal  range  queries.  The  grid  file  [NHS  &4]  Is  a  structure  for  storing  multidimensional  point  data 
designed  to  allow  the  evaluation  of  irregularly  shaped  search  regions  in  such  a  way  that  the  complexity  of 
the  region  affects  CPU  time  but  not  disk  accesses.  The  latter  limit  the  performance  of  a  data  base 
Implementation. 

The  grid  file  partitions  space  Into  raster  cells  and  assigns  data  buckets  to  cells.  The  partition  information  is 
kept  in  scales,  one  for  each  axis  of  space;  the  assignment  is  recorded  In  an  array  called  grid  directory.  The 
directory  Is  likely  to  be  large  and  must  therefore  be  kept  on  disk,  but  the  scales  are  small  and  can  be  kept  in 
central  memory.  Therefore,  the  grid  file  reallies  the  two-disk-access  principle  for  single  point  retrieval 
(exact  match  query):  by  searching  tJie  scales,  the  k  coordinates  of  a  data  point  are  converted  into  Interval 
indices  without  any  disk  accesses;  these  indices  provide  direct  access  to  the  correct  element  of  the  grid 
directory  on  disk,  where  the  bucket  address  Is  located.  In  a  second  access  the  correct  data  bucket  (l.e.  the 
bucket  that  contains  the  data  point  to  be  searched  for.  If  it  exists)  Is  read  from  disk.  A  query  region  Q  Is 
matched  against  the  scales  and  converted  Into  a  set  I  of  index  tuples  that  refer  to  entries  in  ^e  directory. 
Only  after  this  preprocessing  do  we  access  disk  to  retrieve  the  correct  pages  of  the  directory  snd  the  correct 
data  buckets  whose  regions  Intersects  Q  (Fig.  2.6), 


Fig.  2.6:  Query  in  a  grid  file. 

A  geometric  join  query  Is  answered  In  an  analogous  way.  Let  f  and  f  be  the  two  grid  files  Involved,  and  let 
H  and  //'  be  the  underlying  hlgher-dimeralonal  spaces.  The  scales  of  f  and  P  define  a  grid  on  the  Cartesian 
product  H  X  H'.  The  cells  of  tl.ls  grid  which  Intersect  the  search  region  In  H  X  H'  are  determined  by 
matching  the  scales  of  th.e  two  grid  files  against  the  search  region.  As  in  the  case  of  proximity  queries  on  a 
s’.ngle  grid  file  this  computation  needs  no  access  to  disk.  If  a  cell  Intersects  the  search  region  the 
corresponding  pair  of  buckets  (Bf  ,  Bf)  Is  axxessed  from  disk  via  the  grid  directories  of  f  and  P.  If  the 
Cartesian  product  of  Ih.e  bucket  regions  of  Bf  and  Bp  Is  completely  contained  In  the  search  region  all  pain 
cf  objects  correspondirg  to  pairs  of  points  (pf,  pf)  with  pf  €  Bf  and  pp  €  Bp  fUlflll  the  Join  condition.  If 
tb.e  Cartesian  product  of  the  bucket  regions  of  Bf  and  Bp  Is  not  completely  contained  In  the  search  region 
all  pairs  of  points  (pf  ,  pp)  with  pf  £  Bf  and  pp  £  Bp  must  be  checked  In  order  to  see  whether  they  He 
Inside  or  outside  the  search  region,  1.  e.  whether  the  ccTrespondlng  pain  of  objects  filflll  the  Join  condition. 
A  buffer  of  mlrumal  size  of  two  pages  receives  pain  of  data  buckets  (Bf ,  Bp)  acx  'riing  to  a  scheduling 
policy  similar  to  the  one  rr.enConed  In  (MKY  811. 

2.4  The  arid  fUe  software  package 

Th.e  gild  file  Is  Implemented  In  ModuU-2  and  FORTRAN-77  as  a  portable  data  management  package  of 
a-bout  atrut  58riO  lines  of  source  code.  The  Modula-2  version  (Hln  85a,  85b]  runs  on  the  DEC-VAX  11 
under  VMS,  on  the  DEC-PDP  11  under  RT-11  and  on  some  personal  computen  based  on  the  Motorola 
0-VXX)  processor.  The  FORTRA.N  77  version  has  been  developed  on  a  DEC-VAX  11.  The  package  Includes 
a  Prolog  Interpreter  that  gives  the  user  Interactive  access  to  the  data  store  '  In  grid  files,  and  serves  as  a 
powerful  query  language  that  permits  deduction. 


Fig.  2.4b:  Search  region  for  an  Intersection  query  with  a  line  L. 

Geometric  Join  query.  Let  Q*  be  another  class  of  simple  spatial  objects  with  parameter  space  H’,  and  T  C 
O'.  For  every  A  €  0  let  C  H’be  the  set  of  all  points  In  fT  representing  A’  C  0*  such  that  A  and  A’ 
Intersect  Denote  by  the  point  In  H  representing  a  spatial  object  A  €  0.  The  region  In  the  Cartesian 
product  H  y.  H'  that  contains  all  points  representing  pairs  (A,  A’)  €  F  X  T  of  intersecting  objects  Is  the 
union  of  the  sets  {Py^}  X  H\  for  all  A  €  O;  this  region  Is  particularly  simple  for  the  different  classes  of 
simple  spatial  objects. 

Let  Q  be  the  class  of  points.  O’  the  class  of  Intervals  on  a  straight  line.  Then  //  X  Is  the  3-dlmenslonal 
space.  All  pairs  (p,  I)  of  points  p  with  coordinate  x  and  Intervals  I  =  (cx,  dx)  such  that  p  €  I  are  represented 
by  points  lying  in  the  solid  shown  In  Fig.  IS.  This  solid  Is  obtained  by  moving  the  search  region  for  a 
polnt-ln-lnterval  query  along  the  bisector  in  the  x-cx-plane. 


dx 


search  region  for 
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Fig.  2S:  Search  region  for  a  geometric  Join  query. 
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Fig.  2.2:  Search  region  for  a  point  query  In  the  class  of  aligned  rectangles  In  the  plane. 

3)  Let  0  be  the  class  of  circles  In  the  plane.  As  parameters  for  the  representaaon  of  a  circle  as  a  point  In 
3-cllmensJonal  space  we  choose  the  coordinates  of  Its  center  (cx,  cy)  and  Its  radius  r.  All  circles  which 
overlap  a  point  q  are  represented  In  the  corresponding  3-dlmensloral  space  by  points  lying  In  the  cone  with 
vertex  q  shown  In  Fig,  2,4a.  The  axis  of  the  cone  Is  parallel  to  the  r-axls  (the  extension  parameterX  Its 
vertex  Is  q  considered  as  a  point  In  the  cx-cy-plane  (the  subspace  of  the  location  parameters). 

Point  set  query.  Given  a  set  Q  of  points,  the  region  In  //  that  contains  all  pwlnts  representing  objects  A  €  F 
which  Intersect  Q  Is  theunJon  of  the  reglorj  in  H  that  result  from  the  point  queries  for  each  point  in  Q.  The 
urJon  o/cones Is  a  particularly  simple  region  In  //If  the  query  set  Q  Is  a  simple  spatial  object. 

1)  Let  n  be  the  class  of  Intervals  on  a  straight  line.  An  Interval  I  =  (cx.  dx)  Intersects  a  query  Interval  Q  = 
(cq.  dq)  If  and  only  If  Its  representing  point  lies  In  the  shaded  region  shown  In  Fig.  2J:  this  region  Is  given 
by  the  inequalities  cx  -  dx  <  cq  +  dq  and  cx  +  dx  ^  cq  -  dq. 

2)  Let  Q  be  the  class  of  aligned  rectangles  In  the  plane.  If  Q  Is  also  a.n  aligned  rectangle  then  Q  Is  again 
treated  as  the  Cartesian  product  of  two  classes  of  Intervals,  one  along  Lhe  x-axls,  the  other  along  the  y-axls. 
All  rectangles  which  Intersect  Q  are  represented  by  points  In  4-dimenslonaJ  space  lying  In  the  Cartesian 
product  of  two  Interval  Intersection  query  regions. 

3)  Let  0  be  the  class  of  circles  in  the  plane.  All  circles  which  Intersect  a  line  segment  L  are  represented  by 
points  lying  In  the  cone-shaped  solid  shown  In  Fig.  2.4b.  This  solid  Is  obtained  by  embedding  L  In  the 
cx-cy-plane.  the  subspace  of  the  location  parameters,  and  moving  the  cone  with  vertex  at  q  along  L. 


Intervals  on  a  long  line  clustering  along  the  diagonal,  leaving  large  regions  of  a  large  en^beddlng  space 
unpopulated:  whereas  the  same  set  of  Intervals  represented  by  a  location  parameter  cx  and  an  extension 
parameter  dx.  fills  a  smaller  embedding  space  In  a  much  more  uniform  way.  With  the  assumption  of 
bounded  d.  this  data  distribution  Is  easier  to  handle. 
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11  Intersection  queries  lead  to  cone-shaped  search  regions 

InterswUon  Is  a  basic  component  of  other  proximity  queries,  and  thus  deserves  special  attenUen.  C/JD 
design  rules,  for  example,  often  require  different  objects  to  be  separated  by  some  minimal  distance.  This  is 
equivalent  to  requiring  that  objects  surrounded  by  a  rim  do  not  intersect.  Given  a  class  0  of  simple  spatial 
objects  wi’di  parameter  space  //,  and  a  set  T  C  Q  of  simple  objects  represented  as  points  In  H,  we  consider 
three  Types  of  queries; 

-  point  query:  given  a  query  point  q,  find  all  objects  A  €  T  for  which  q  £  A. 

-  point  set  query;  given  a  set  Q  of  points,  find  all  objects  ACT  which  Intersect  Q. 

-  geometric  Join  query:  given  arolher  class  O'  of  spatial  objects  with  parameter  space  H\  and  a  set  P  C  O', 

find  all  pairs  (A.  A’)  €  F  X  F  of  Intersecting  objects. 

Point  query.  For  a  query  point  q  compute  the  region  In  H  that  contains  all  points  representing  objects  In  F 
which  overlap  q. 

1)  Let  Q  be  the  class  of  Intervals  on  a  straight  line.  An  Interval  given  by  Its  center  cx  and  Its  half  lcr.gth  dx 
overlaps  a  point  q  with  coordinate  qx  If  and  only  If  cx  -  dx  <  qx  ^  cx  +  dx. 

2)  The  class  0  cf  aligned  rectungles  In  the  plar,e  (with  parameters  cx,  cy,  dx,  dy)  can  be  treated  as  the 
Curiesian  product  of  two  classes  cf  Intervals,  one  along  the  x-axls,  the  other  along  the  y-axls.  All  rectangles 
which  contain  a  given  point  q  are  represented  by  points  In  4 -dimensional  space  lying  in  the  Cartesian 
product  of  two  polnt-ln-lntcrval  query  regions  (Fig.  2.2).  The  region  Is  shown  by  Its  projections  Into  the 
cx-dx-plane  and  the  cy-dy-plane. 


simple  object  such  as  a  point  or  a  line  segment  Even  If  the  solid  and  the  query  are  far  apart,  all  the 
components  of  the  solid  must  be  examined  in  a  tree  traversal  to  detect  this.  What  Is  lacking  Is  some  concisely 
stat^  geometric  Information  that  describes  global  properties  of  the  solid  and  Its  location  In  space. 

2  An  approach  to  combined  geometric  modeling  and  storing: 

Approximation,  transformation  'n  parameter  space,  grid  file 

The  technique  we  now  present  for  modeling  and  storing  spatial  objects  la  based  on  1)  approximation  of 
complex  spatial  objects  by  simple  shapes,  eg.  contalnerx  2)  transformation  of  simple  spatial  objects  into 
points  in  higher-dlmensioruil  parameter  spaces  and  3)  the  grid  file  for  point  Storage. 

Complex,  Irregularly  shaped  spatial  objects  can  be  represented  or  approximated  by  simpler  ones  In  a 
variety  of  ways,  for  example:  decomposition,  as  In  a  quad  tree  tessellation  of  a  figure  Into  dl^olnt  raster 
squares  of  site  as  large  as  possible;  representation  as  a  cover  of  overlapping  simple  shapes;  enclosing  it  in  a 
container  chosen  from  a  class  of  simple  shapes.  The  container  technique  allows  efficient  processing  of 
proximity  queries  because  It  preserves  the  most  important  properties  for  proxlmlty-based  access  to  spatial 
objects.  In  particular  It  does  not  break  up  the  object  Into  components  that  must  be  processed  separately, 
and  It  eliminates  many  pwtenllal  tests  quickly  (If  two  containers  don’t  Intersect,  the  objects  within  won’t 
either).  As  an  example,  consider  finding  all  polygons  that  Intersect  a  given  query  polygon,  given  that  each 
of  them  Is  enclosed  In  a  simple  container  such  as  a  circle  or  an  aligned  rectangle.  Testing  two  polygons  for 
intersection  Is  an  expensive  operation  as  compared  to  testing  their  containers  for  Intersection.  The  cheap 
container  test  excludes  most  of  the  polygons  from  an  expensive,  detailed  Intersection  check. 

Any  approximation  techrJque  limits  the  primitive  shapes  that  must  be  stored  to  one  or  a  few  types  for 
example  aligned  rectangles  or  boxes.  An  Instance  of  such  a  type  Is  determined  by  a  few  parameters,  such  as 
coordlrates  of  its  center  and  its  extension,  and  can  be  considered  to  be  a  point  In  a  (higher -dimensional) 
parameter  spuce.  This  tmrs formation  reduces  object  storage  to  point  storage.  Increasing  the  dimensionality 
of  the  problem  but  without  loss  cf  information.  Combined  with  an  efficient  multidimensional  data  structure 
for  point  storage  it  Is  the  basis  for  an  effective  Implementation  of  data  bases  of  spatial  objects, 

2.1  Transformation  to  parameter  space 

Cor.slder  a  class  of  simple  spatial  objects,  such  as  aligned  rectangles  In  the  plane  (i.e.  with  sides  parallel  to 
the  axes).  Within  Its  class,  each  object  is  defined  by  a  small  number  of  parameters.  For  example,  an  aligned 
rectar.gle  is  determined  by  its  center  (cx,  cy)  and  the  half-length  of  each  side,  dx  and  dy. 

An  object  defined  within  Its  class  0  by  k  parameters,  can  be  considered  to  be  a  point  In  a  k-dlmenslonal 
parameter  space  H  assigned  to  Q.  For  example,  an  aligned  rectangle  becomes  a  point  In  4-dlmenslorral 
space.  All  of  the  geometric  and  topological  properties  of  an  object  can  be  deduced  from  the  class  It  belongs  to 
and  from  the  coordinates  of  its  corre.^pondtng  point  In  parameter  space. 

Different  choices  of  the  parameter  space  H  for  the  same  class  0  of  objects  are  appropriate,  depending  on 
characteristics  of  the  data  to  be  processed.  Some  considerations  that  may  determine  the  choice  of 
parameters  axe: 

1)  Distinction  between  location  parameters  and  extension  pammeters.  For  some  classes  of  slm.ple  objects  It 
Is  reasonable  to  distinguish  location  parameters,  such  as  the  center  (cx,  cy)  of  an  aligned  rectangle,  from 
extension  pwrameters,  such  as  the  half  sides  dx  and  dy.  This  distinction  Is  always  possible  for  objects  that 
can  be  described  as  Cartesian  products  of  sphere  of  various  dimensions.  For  example,  a  rectangle  Is  the 
product  of  two  1-dimerrsloaaJ  spheres,  a  cylinder  the  product  of  a  1-dlmenslonal  and  a  2-dlmenslonal 
sphere.  Whenever  this  distinction  can  be  made,  cone-shaped  search  regions  generated  by  proximity  queries 
as  described  in  section  2-3  have  a  simple  Intuitive  Interpretation:  The  subspace  of  the  location  parameters 
acts  as  a  "mirror'’  that  reflects  a  query. 

2)  Independence  of  parameters  uniform  distribution.  As  an  example,  corsidcr  the  class  of  all  Intervals  on  a 
straight  line  (Fig.  2.1).  If  Intervals  are  represented  by  their  left  and  right  endpoints,  lx  and  rx,  the  constraint 
lx  <  rx  restricts  all  representations  of  these  Intervals  by  points  (lx,  rx)  to  the  triangle  above  the  diagonal. 
Any  data  structure  that  orgarUzes  the  embedding  space  of  the  data  points,  as  opposed  to  the  particular  set  of 
points  that  must  be  stored,  will  pay  some  overhead  for  representing  the  unpopulated  half  of  the  embedding 
space.  A  coordinate  transformation  that  distributes  data  all  over  the  embedding  space  leads  to  more 
efficient  storage.  The  phenomenon  of  nonunlform  data  distribution  can  be  worse  than  this,  in  most 
appllcallorvs,  the  building  blocks  from  which  complex  objects  are  built  are  much  smaller  than  the  space  in 
w  hich  they  arc  embedded,  as  the  size  of  a  brick  Is  small  compared  to  the  size  of  a  house.  If  so.  parameters 
such  as  lx.  rx  that  locate  boundaries  of  an  object,  arc  highly  dependent  on  each  other.  Fig.  2.1  shows  short 


1 J  The  conTcntlonal  data  base  approach  to  "aoa'Staodard**  data 

Data  base  technology  has  developed  over  the  past  two  decades  In  response  to  the  needs  of  commerdal 
data  processing.  The  key  concepts  Introduced  and  sup>ported  by  data  base  software  mirror  the  reality  that 
used  to  be  handled  manually  by  ofTlce  clerks.  Large  quantities  of  records  of  a  few  dlfTerent  types. 
Identified  by  a  small  number  of  attributes,  mostly  retrieved  In  response  to  relatively  simple  queries:  point 
queries  that  ask  for  the  presence  or  absence  of  one  particular  record,  interval  or  range  queries  that  ask  for 
all  records  whose  attribute  values  lie  within  given  tower  and  upper  bounds.  More  complex  queries  tend  to 
be  reduced  to  these  basic  typ^ 

Data  base  software  has  yet  to  take  Into  account  the  specific  requirements  of  geometric  computation,  as  can 
be  seen  from  the  termlrralogy  used:  Geometric  objects  are  lumped  Into  the  amorphous  pool  of 
"non-standard"  applications.  The  sharp  distinction  between  the  logical  view  presented  to  the  user  and  the 
physical  aspects  that  the  Implementer  sees  has  been  possible  In  conventional  data  base  applications 
because  data  structures  that  allow  efTldent  handling  of  point  sets  are  well  understood.  The  same  dlstlction 
Is  premature  for  geometric  data  bases:  In  interactive  applications  such  as  CAD  cfndency  Is  the  real  Issue, 
and  until  we  understand  geometric  storage  techrvlques  better  we  may  not  be  able  to  afford  the  luxury  of 
studying  geometric  modeling  divorced  from  physical  storage.  Consider  the  following  example. 

A  set  of  polyhedra  might  be  stored  In  a  relational  data  base  by  using  the  bouixiary  representation  (BR) 
approach:  a  polyhedron  p  Is  given  by  Its  faces,  a  face  f  by  its  bounding  edges,  an  edge  e  by  its  endpoints  sj 
and  %2-  Four  relations  polyhedra,  faces,  edges  and  po/nXs  might  have  the  following  structure: 

A  tuple  In  the  relation  polyhedra  Is  a  pair  (pj.  f^J  of  Identifiers  for  a  polyhedron  and  a  face:  fjj  is  a  face  of 
polyhedron  pj. 

A  tuple  In  the  relation  faces  Is  a  pair  (fi{,  ej)  of  Identifiers  for  a  face  and  an  edge:  ^  Is  a  bounding  edge  of 
face  fjj. 

A  tuple  In  the  relation  edges  Is  a  triple  (^.  sj,  of  Identifiers  for  one  edge  and  two  points;  S]  and  s^,  are 
the  endpoints  of  edge  e^. 

A  tuple  in  the  relation  po/n/r  is  a  triple  (s^.  x.  y):  s^  Is  the  Identifier  of  a  point,  x  and  y  are  Its  coordinates. 

This  representation  smashes  an  object  Into  parts  which  are  spread  over  different  relations  and  therefore 
over  the  storage  medium.  The  question  whether  a  polyhedron  P  Intersects  a  given  line  1  Is  answered  by 
IntersecUng  each  face  f^  of  a  polyhedron  pj  with  1.  If  the  tuple  (pj,  fi^)  In  the  relation  pelyhedra  contains 
the  equsUon  of  the  corresponding  plane,  the  Intersection  point  of  the  plane  and  the  line  1  can  be  computed 
without  accessing  other  relations.  But  In  order  to  determine  whether  this  Intersection  point  lies  Inside  or 
outside  the  face  fv  requires  accessing  tuples  of  edges  and  polnls,  1.  e.  accessing  different  blocks  of  storag^ 
resulUng  In  many  more  disk  accesses  than  the  geometric  problem  requires. 

1.4  Geometric  modeling  separated  from  storage  considerations 

In  this  early  stage  of  development  of  geometric  data  base  technology,  we  cannot  afford  to  focus  on 
modellrg  to  the  exclusion  of  implementation  aspects.  In  graphics  and  CAD  the  red  Issue  Is  efp.ctency. 
1/10-th  of  a  second  Is  limit  of  human  time  resolution,  and  a  designer  works  at  maximal  efndency  when 
"trivial"  requests  are  displayed  "Ir^tanlaneously".  This  allows  a  couple  of  disk  accesses  only,  which  means 
that  geometric  ard  other  spatial  attributes  must  be  part  of  the  retrieval  mechanism  If  common  geometric 
q'ueries  (Inlersectlcn.  Inclusion,  point  queries)  are  to  be  handled  effldently. 

A  key  problem  that  affects  efP.dency  Is  how  to  reduce  complex  objects  to  simpler  ones  chosen  from 
predefined  primitives.  Among  the  standard  techniques  known  we  have  already  discussed  how  boundary 
representations  stored  in  a  relational  data  base  prevent  effident  access  based  on  geometric  queries.  The 
problem  of  an  object  being  tom  apart  happens  also  In  another  standard  modeling  technique,  constructive 
solid  geometry  CSG.  Let  us  briefly  discuss  the  consequences  of  basing  the  physical  storage  structure 
directly  on  such  modeling  techniques. 

In  constructive  solid  modeling  a  complex  object  Is  constructed  from  simple  primlUves,  such  as  cubes  or 
spheres,  by  means  of  Boolean  operations  union,  intersection  and  dlfTerencc.  The  construction  process  Is 
represented  by  a  tee.  Each  leaf  of  a  CSG  tree  contains  a  simple  object,  each  Internal  node  contains  a 
Boolean  operation.  To  each  node  a  geometlc  tansformatlon  such  as  scaling,  translation  and  relation  may 
be  assigned.  The  Boolean  operation  Is  applied  to  the  objects  represented  by  the  left  and  right  subtree  of 
the  node.  A  geometric  ta.ns formation  assigned  to  a  leaf  Is  applied  to  the  simple  object  stoted  In  the  leaf,  a 
gecm.e'jic  transformation  assigned  to  an  Internal  node  Is  applied  to  the  object  resulting  from  the  Boole.m 
operation  stored  In  this  node.  Now  consider  the  query  whether  a  solid  In  CSG  representaUon  Intersects  a 
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The  first  generation,  characterized  by  numerical  ccmpuling,  led  to  the  development  of  many  new 
algorithms.  It  transformed  numerical  analysis  from  a  craft  to  be  practiced  by  every  applied  mathematician 
Into  a  field  for  specialists.  It  soon  became  obvious  that  writing  good  (efRclent.  robust)  numerical  software 
requires  so  much  knowledge  and  effort  that  this  task  cannot  be  left  to  the  applications  programmer.  The 
development  of  large  portable  numerical  libraries  became  one  of  the  major  tasks  for  professional 
numerical  analysts. 

The  second  generation,  hatched  by  the  needs  of  comirerdal  data  processing,  led  to  the  development  of 
many  new  data  structures.  It  focused  attention  on  the  problem  of  efficient  management  of  large,  dynamic 
data  collections.  Initially  under  batch  processing  conditions.  Searching  and  sorting  were  recognized  as 
basic  operations  whose  time  requirements  turned  out  to  be  the  bottleneck  for  many  applications.  Data 
base  technology  emerged  to  shield  the  end-user  from  the  details  of  Implementation  (storage  techniques, 
features  dependent  on  hardware-  and  operating  systemX  by  presenting  the  data  In  the  form  of  logical 
models  that  highlight  relationships  among  data  Items  rather  than  their  Internal  representation,  and  by 
Introducing  the  abstraction  of  access  path  to  hide  detailed  access  algorithms  of  underlying  data  structures. 

We  are  now  on  the  threshold  of  a  third  generation  of  applications,  dominated  by  computing  with  pictorial 
and  geometric  objects.  This  change  of  emphasis  Is  triggered  by  toda/s  ubiquitous  Interactive  use  of 
personal  computers,  and  their  increasing  graphics  capabilities.  It  Is  a  simple  fact  that  people  absorb 
information  fcistest  when  it  Is  presented  In  pictorial  form,  hence  computer  graphics  and  the  underlying 
processing  of  geometric  objects  will  play  a  role  in  the  majority  of  computer  applications.  The  Held  of 
computational  geometry  has  emerged  as  a  scientific  discipline  during  this  past  decade  In  response  to  the 
growing  importance  of  processing  pictorial  and  geometric  objects.  It  has  already  created  novel  and 
Interesting  algorithms  and  data  structures,  and  Is  beginning  to  Impact  data  base  technology  under  the  label 
(hopefully  temjxjrary)  of  non-standard  database  applications.  In  order  to  understand  how  geometric 
computation  is  likely  to  affect  data  bases.  It  Is  useful  to  survey  some  milestones  In  this  rapidly  developing 
field  which  Is  replacing  the  traditional  areas  of  numeric  computation  and  of  data  management  as  the  major 
research  topic  in  algorithm  analysis. 

1.2  Computational  geometry  *  theory  and  practice 

During  the  seventies  geometric  problems  caught  the  attention  of  researchers  In  concrete  complexity 
theory.  They  brought  to  bear  the  finely  honed  tools  of  algorithm  analysis  and  achieved  rapid  progress. 
Bementary  problems  (e.  g.  determining  intersections  of  simple  objects  such  as  line  segments,  aligned 
rectangles,  polygons)  yielded  elegantly  to  general  algorithmic  principles  such  as  divide- and-conquer  or 
plane  sweep.  But  In  many  Instances  a  surprisingly  large  Increase  of  difficulty  showed  up  In  going  from  bwo 
to  three  dimensions:  for  example.  Intersection  of  polyhedra  Is  still  an  active  research  topic  where  major 
efficiency  gains  are  to  be  expected.  The  theory  of  computational  geometry,  although  well  underway,  has  as 
yet  explored  only  a  fraction  of  Its  potential  territory. 

The  practice  of  computational  geometry  Is  even  less  well  understood.  Many  Important  geometric  problems 
in  computer-aided  dcsiga  In  geographical  data  processing.  In  graphics  do  not  lend  themselves  to  being 
studied  and  evaluated  by  the  asymptotic  perfotmance  formulas  that  the  algorithm  analyst  cherishes.  For 
exa-mple,  asymptotics  does  not  help  In  answering  the  question  whether  we  can  access  an  object  In  one  disk 
access  or  two,  thus  being  able  to  display  It  "Instantaneously"  on  the  designer's  screen  -  realistic 
as3un-ipaor,s  about  the  size  of  today's  central  memories  are  needed  Nor  will  asymptotics  settle  the 
argument  raging  In  th.e  CAD  community  between  proponents  of  boundary  representaUor.s  and  adherents 
of  constructive  solid  geometry  -  taste,  experience,  and  type  of  application  are  the  relevant  parameters. 
And  below  the  highly  visible  issues  of  object  representa’Jon,  data  structures  and  algorllhm.s  hide  the 
tantalizing  details  of  the  numerics  of  computaUonal  geometry,  such  as  the  problems  caused  by  "braiding 
straight  lines",  which  may  intersect  repeatedly. 

Commercially  available  software  In  computer  graphics  and  CAD  has  not  yet  taken  Into  account  the  results 
of  ccmputatJoruil  geometry.  Stralghtfor-vard  algorithm, s  are  being  used  whose  theoretical  efficiency  Is  poor 
as  compared  to  known  results.  Perhaps  the  straightforward  algorithms  are  better  In  practice  than 
theoretically  optlrruil  ones,  but  such  difficult  questions  have  hardly  been  Investigated,  as  CAD  systems 
development  today  Is  so  labor  Intensive  that  all  resources  are  absorbed  by  just  getting  the  system  to  work, 
and  algorithm  analysis  has  so  far  largely  restricted  Itself  to  theoretically  measurable  performance. 

We  know  by  analogy  with  numerical  anal’, -sis  what  the  next  step  should  be  In  the  maturing  process  of 
cxampulatlonal  geometry:  The  developmenl  of  efficient,  portable,  robust  program  libraries  for  the  most  basic, 
frequent  geometric  operations  on  stardard  rcprcorntallons  of  geometric  objects.  In  other  words,  we  must 
develop  the  geometric  subroutine  library  of  CAD.  thus  exposing  theoretl^  results  to  stringent  practical 
tests. 
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Fig.  2;  Conceptual  Scheme  of  the  CSO-Approach. 


In  the  classical  relational  model,  data  is  organized  as  record  instances  (tuples)  in  tables 
(relations).  There  are  no  schema  defined  relationships  such  as  PART-OF  and  IS-A 
structures.  Instead,  a  high  level  language  (predicate  calculus  or  relational  algebra)  is  used  for 
exploiting  relationships  based  on  values.  However,  the  mapping  of  highly  interrelated  data 
into  tuples  in  one  or  more  relations  has  to  be  done  entirely  by  the  user.  For  instance,  the 
important  feature  of  a  generic  structure  which  says  that  all  descendant  objects  must  bear  the 
same  key  domains  as  their  ascendants  has  to  be  enforced  by  the  user  himself.  The  relational 
model  does  not  allow  individual  objects  to  be  uniformly  referred  to  regardless  of  the  generic 
class  in  which  they  appear. 

Other  drawbacks  using  the  classical  relational  model  are  due  to  normalization.  A  relation  is 
said  to  be  in  First  Normal  Form  [UUman  1982]  if  and  only  if  it  satisfies  the  constraint  that  it 
contains  atomic  values  only.  Note  that  in  our  example  of  a  CSG-scheme,  this  condition  is 
very  inconvenient  For  instance,  a  TRANSFORMED-PART  may  be  described  by  a  part 
number  and  a  transformation  matrix,  e.g.  a  4x4-matrix  in  homogeneous  coordinates.  To 
describe  these  facts  in  a  normalized  relation,  sixteen  attributes  must  be  introduced 
artificially  which  corre.spond  to  the  matrix  arguments.  Of  course,  one  would  instead  only 
define  three  arguments  PARAMl,  PARAM2  and  PARAM3  for  motion  parameters,  plus  a 
further  attribute  MOTION  which  denotes  if  it  is  a  translation,  a  rotation,  or  a  scaling 
operation.  In  any  case,  the  First  Normal  Form  is  cumbersome. 

In  conclusion,  the  study  of  the  CSG-approach  in  solid  modeling  suggests  the  following 
extensions  to  the  relational  model:  There  should  be  a  way  of  defining  PART-OF  and  IS-A 
structures  e.xplicitly  to  the  system  in  order  to  give  the  user  the  possibility  of  querying  an 
object  or  pan  of  it  as  a  whole  rather  than  assembling  different  relations  and  thinking  about 
all  ’'nown  interrelationships.  In  addition,  the  First  Normal  Form  should  be  dropped.  Or  at 
least,  the  user  should  have  a  direct  way  of  storing  matrices  as  data  types  into  a  tuple,  and  the 
database  system  should  incorporate  features  for  non-atomic  fields  into  its  calculus  or 
algebra. 


22.  Boundary  Representation 

With  the  Boundary  Representation  model,  solids  are  described  by  a  collection  of  faces  which 
in  turn  are  represented  by  their  bounding  edges  and  faces.  So  called  Euler  operators  allow 
incremental  manipulation  of  the  objects  while  restoring  well-formedness  of  the  surfaces: 
closedness,  orientability,  non.self-intersection.  boundedness,  and  connectivity  [Eastman  and 
Weiler  1979].  These  topological  properties  are  condensed  in  the  Euler- Poincar6  formula: 
with  f  faces,  e  edges,  v  vertices,  r  inner  loops  of  faces  called  rings,  c  cavities  or  hollow  tubes, 
and  g  holes  through  the  body  (or  genus  g  corresponding  to  the  number  of  handles  in  graph 
topology)  the  following  condition  holds: 

f-e-*-v-r  =  2*(c-g) 

The  practical  relevance  of  the  formula  is  emsuring  that  shapes  are  topologically  well-formed; 
e.g.  its  application  eliminates  the  danger  of  ill-formed  solids  such  as  the  Klein  bottle.  If  we 
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consider  the  above  formula  as  a  hyperplane  In  six-dimensional  space,  the  law  restricts  the 
valid  transitions  to  a  subset  of  all  those  comblnatoiically  possible.  Of  course,  the  desired  set 
of  Euler  operators  should  cover  the  hyperplane;  a  possible  spanning  set  of  five  primitive 
operators  may  be  defined  as  follows: 


MEF  resp.  KEF 
MEV  resp.KEV 
MEKR  resp.  KEMR 
MFVC  resp.  KFVC 
MFKGR  resp.  KFMGR 


f  e  V  r  c  g 
1  1  0  0  0  0 
0  110  0  0 
0  1  0-1  0  0 
10  10  10 
1  0  0-1  0-1 


The  Euler  operator  MEF  stands  for  ’’Make  Edge  and  Face”  which  obviously  does  not  change 
the  above  characteristic,  it  is  also  invers  to  KEV,  i.e.  "Kill  Edge  and  Face”.  Any  transition  in 
the  Eulerian  plane  can  now  be  represented  as  a  linear  combination  of  the  five  primitive 
Euler  operators.  Each  of  these  or  a  combination  enables  the  construction  of  a  possible 
unique  topology.  They  reduce  bookkeeping  requirements  needed  to  guarantee  that  the 
resulting  shapes  are  well-formed,  i.e.  non-intersecting,  closed,  and  orientable. 

We  now  discuss  the  conceptual  scheme  of  solids  described  in  boundary  representation  (see 
Fig.  3).  The  structure  of  a  bounded  shape  model,  i.e.  OBJECT,  is  comprised  of  spatial 
surfaces  named  faces.  Each  FACE  is  bounded  by  one  or  more  loops  of  edges  where  each 
loop  is  the  concatenation  of  line  segments,  Le.  edges,  into  a  closed  RING.  EDGES  are 
bounded  by  VERTICES  at  their  intersections;  in  our  scheme,  every  edge  is  given  by  a 
START  and  END  vertex,  and  it  topologically  points  to  the  LEIT  and  RIGHT  ring 
respectively. 


Fig.  3:  Conceptual  Scheme  of  the  BR-Approach. 


Using  a  DBMS  for  storing  objects  in  boundary  representation  is  advantogeous  for  the 
following  reasons.  For  instance,  the  discussed  Euler  operators  are  atomic  in  the  sense  that 
they  topologically  guarantee  to  handle  manifolds  consistently,  i.e.  to  fulfill  the 
Euler- Poincar6  formula.  To  construct  the  topology  of  a  cube  for  instance,  a  sequence  of  a 
MFVC,  seven  MEV,  and  five  MEF  operators  is  needed.  Thus,  the  notion  of  transaction  in 
database  theory  [Ullman  1982]  might  be  very  helpfui  to  better  control  coasistency:  The 
sequence  of  Euler  operators  would  start  with  a  BEGIN  TRANSACTION  command  and 
finish  with  an  END  TRANSACTION  command.  The  transaction  mechanism  is  such  that 
other  transactions  (or  users)  do  not  see  the  changes  of  a  transaction  until  this  transaction 
commits.  When  it  commits,  the  whole  sequence  of  Euler  operators,  e.g.  the  topology  of  a 
cube,  become  visible  to  other  users.  If  a  transaction  does  not  commit  but  aborts  or 
terminates  abnormally,  any  change  to  the  data  are  undone  and  other  transactions  will  see 
none  of  the  changes. 

There  are  some  drawbacks  when  using  relational  database  technology  for  bounded  surfaces. 
First,  the  user  is  forced  to  define  keys  such  asO^,F<!^.R/!f,EXf  and  V;!^.  These  values  are 
necessary  to  uniquely  identify  each  tuple  within  a  relation.  However,  to  enumerate  all  faces 
of  an  object,  all  rings  of  a  face,  all  edges  of  a  ring  or  both  vertices  of  an  edge  should  be 
superfluous  when  interacting  through  the  graphical  interface  of  the  solid  modeler. 


Each  instance  of  the  BR-scheme  is  a  group  of  tuples  comprising  a  single  root  tuple  which 
defines  the  object,  and  several  dependent  tuples  in  dikinct  relations  which  form  its 
boundary,  l.e.  faces,  rings,  edges,  and  vertices.  Even  if  a  structure  mch  as  PART-OF  would  be 
expressed  relationally  in  terms  of  matching  values,  it  could  not  be  manipulated  as  a  single 
object  For  example,  in  order  to  delete  an  object  the  user  must  issue  one  delete  statement 
for  each  tuple  in  all  dependent  relations  of  the  object 

Also  the  First  Normal  Form  condition  brings  disadvantages  for  a  BR-scheme.  For  instance,  a 
vertex  may  be  described  by  its  positional  number  and  coordinates.  Due  to  normalization,  the 
coordinates  may  not  be  treated  as  vectors  but  have  to  be  distributed  into  three  attribute 
domains,  namely  for  x-,  y-,  and  z-values. 


3.  A  SURROGATE  MODEL 


We  now  describe  a  surrogate  concept  as  the  basis  for  an  engineering  database  system,  we 
demonstrate  the  usability  of  surrogates  for  defining  PART-OF  and  IS-A  structures  in  solid 
modeling,  and  we  give  some  data  retrieval  and  manipulation  considerations. 


3.1.  Sonogates  versus  User  Keys 

In  [Hall  et  al.  1976],  it  was  pointed  out  that  the  relational  model  cannot  denote  an  individual 
object  independently  of  its  attributes.  In  other  words,  what  would  happen  if  a  particular  part 
number  (unique  identifier)  in  a  CAD  database  is  replaced  by  a  new  one;  Is  it  a  change  to 
that  part  number  or  a  replacement  by  a  new  part  with  the  same  characteristics?  To  solve  this 
problem.  Hall  et  al.  propose  to  use  surrogates  as  "data  model  representatives"  of  the  entities 
I  •  (unlike  tuple  identifiers  used.  e.g.  in  System  R  [Astrahan  et  aL  1976])  and  draw  the  following 
distinction: 

SURROGATES  act  as  invariant  values  for  individual  entities;  these  values  can  appear 
at  different  places  in  the  database  to  link  entities  together. 

USER  KEYS  act  as  unique  identifiers  under  user  control  to  identify  an  individual 
entity. 

One  extension  of  the  relational  model  [Codd  1979]  suggests  a  unary  relation  for  each  entity 
type  to  list  all  the  surrogates  of  entities  which  are  currently  recorded  in  the  database.  Codd’s 
entity  integrity  constraint  allows  insertions  and  deletions  of  surrogate  values  but  not  updates 
and  null  values. 

Surrogates  can  be  used  to  provide  both  fast  access  and  storage  independence.  Deen's 
implementation  [Deen  1982]  employs  key  compression  to  provide  a  more  uniform 
distribution,  a  hashing  algorithm  to  place  tuples  on  data  pages,  and  an  indexing  technique  to 
allow  fast  sequential  access.  However,  Deen’s  surrogates  are  similar  to  tuple  identifiers  and 
are  generated  from  primaiy  keys.  Therefore,  whenever  a  primary  key  value  changes,  its 
corresponding  surrogate  also  changes. 

We  introduce  a  system-controlled  attribute  SURROGATE  and  restrict  the  surrogate  values 
according  to  the  following  rules: 


-  Each  SURROGATE  value  is  system- wide  unique[p.g.  concatenation  of  processor  number, 
database  identification,  and  clock  time  or  sequence  number  per  relation)  in  order  to 
allow  for  merging  of  databases  from  different  sites. 

-  The  values  of  a  SURROGATE  attribute  cannot  be  changed  The  user  has  no  control  over 
the  SURROGATE  values  although  they  may  or  may  not  be  made  available  to  him  (e.g.  it 
seems  appropriate  to  give  surrogate  values  back  to  programmers  as  a  program  variable). 


A  SURROGATE  acts  as  an  invariant  value  for  each  tuple,  and  no  special  attribute  needs  to 
be  chosen  as  the  primary  key.  In  our  example  of  the  BR-approach  for  instance,  the  user 
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could  defme  O#,  F#,  R#,  E#  and  as  surrogates.  The  system  would  then 
automatically  generate  unique  identifiers  whenever  a  tuple  is  inserted  in  a  relation.  In 
addition,  these  values  could  be  used  to  define  the  struaural  semantics  of  the  objects. 

Very  often,  the  user  likes  to  deal  with  user-defined  primary  keys  which  have  some  semantic 
meaning  to  him.  To  avoid  the  introduction  of  two  independent  identifier  concepts,  we 
introduce  a  binding  mecharusm  between  SURROGATES  and  USER  KEYS  via  a  special  index 
called  KEY-INDEX.  This  index  is  restricted  to  a  single  attribute,  i.e.  unique  key,  and  implies 
a  binding  to  its  corresponding  SURROGATE  attribute.  It  is  important  to  note  that  a  user 
key  may  or  may  not  exist  and  may  sometimes  be  changed:  Supporting  access  to  an 
individual  tuple  of  a  relation  is  always  guaranteed  via  the  SURROGATE  values. 

Furthermore,  we  define  two  built-in  ftinctions  to  map  system-generated  SURROGATES 
onto  user  defined  USER  KEYS  and  vice  versa:  KEY(surrogate)  retrieves  the  user  key 
corresponding  to  a  surrogate  value  if  one  exists,  and  SURR(user  key)  retrieves  the  surrogate 
1 1,  value  of  a  specific  user  key.  A  one-to-one  mapping  between  internal  SURROGATE  values 
and  USER  KEYS  is  guaranteed  if  the  attribute  of  the  indexed  column  is  specified  with  a 
NOT  NULL  option.  In  this  case,  both  functions  KEY  and  SURR  yield  a  unique  value  which 
is  never  null  whereas  a  non-existing  operand  produces  an  error  message. 


3J.  PART-OF  and  IS- A  Structures 

Based  on  the  surrogate  concept  introduced  so  far,  we  show  how  the  structural  part  of  our 
conceptual  schemes  for  solid  modeling  can  be  described  more  directly.  In  Fig.  2  for  instance, 
the  entity  set  OBJECT  can  be  referred  to  as  root  relation  which  identifies  its  hierarchical 
subparts.  In  order  to  define  this  hierarchical  structure,  we  introduce  the  new  attribute 
SURROGATE  for  system-generated  values  in  the  root  relation: 

RELATION  Object; 

ATTRIBUTE 
Art#: 

0#: 

Description ; 

IDENT 

Art#  PRIMARY  DOMAIN; 

KEY-INDEX 

Art#.0#; 

END  Object; 

Besides  the  object  number  O^  (as  a  surrogate),  a  user  key  Art/Y  may  be  defined  and 
combined  with  a  KEYTNDEX.  This  index  allows  the  user  to  retrieve  data  by  article 
numbers  rather  than  internal  surrogates.  It  also  m.ay  be  used  to  improve  the  performance  of 
queries  based  on  the  user  key  attribute,  e.g.  when  searching  for  tuples  with  a  given  aiticle 
number. 

The  SURROGATE  columns  have  a  semantic  meaning  besides  technical  properties  such  as 
clustering,  avoiding  composite  keys,  and  improving  performance;  They  may  be  used  to 


Number ; 

SURROGATE; 

StrIngZO; 
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reference  relations.  For  instance,  the  dependent  relation  PART  is  distinguished  by  the 
PART-OF  attribute  that  contains  surrogates  pointing  to  tuples  in  the  corresponding  parent 
relation  OBJECT: 

RELATION  Part; 

ATTRIBUTE 
PM: 

OM: 

Material : 

END  Part; 

Funhermore,  an  additional  column  type  IS-A  may  be  used  to  refer  to  other  relations  which 
correspond  to  a  generalization  hierarchy.  As  an  example,  we  consider  the  relation 
t  •  CYLINDER  which  is  generalized  by  the  relation  PRIMITIVE-PART : 

RELATION  Cylinder; 

ATTRIBUTE 
ZM: 

PM: 

Rad lus ; 

Height: 

IDENT 

CM  PRIMARY  DOMAIN: 

(Radius. Height)  UNIQUE; 

END  Cy 1 1 nder ; 

The  PART-OF  construct  is  used  to  define  hierarchies  of  relations  and  implicitly  expresses  an 
existential  quantifi'  uion;  For  each  instance  of  the  hierarchical  class,  there  exist  objects 
constituting  its  pans.  On  the  other  hand,  the  IS-A  hierarchy  implicitly  expresses  a  universal 
quantification:  Every  instance  of  a  subordinate  class  has  all  the  properties  of  the  more 
general  class.  (Besides  PART  OF  and  IS-A  structures,  an  additional  attribute  type 
RFEERENC’E  OF  may  be  defined  to  refer  to  tuples  of  the  same  or  a  different  hierarchy. 
This  construct,  however,  would  ask  for  specific  semantics,  and  performance  enhancements 
would  become  more  difficult,  i.e.  natural  clustering  of  data  may  no  longer  be  applicable). 


3  J.  Data  Retrieval  and  Manipulations 

To  retrieve  data  from  PART-OF  and  IS-A  hierarchies,  a  user  would  often  have  to  join 
component  relations  with  parent  or  ancestor  relations  which  requires  knowledge  of  the 
external  structure  of  a  complex  object  Instead  of  defining  several  join  predicates  along 
particular  branches  involving  SURROGATE  PART-OF,  and  IS-A  columns,  the  user  may 
specify  an  implicit  join  operator.  By  this,  the  whole  implicit  structure  of  aggregation  or 
generalization  concepts  become  more  transparent  and  may  be  easier  handled  at  the  user 
interface. 


Number ; 

IS-A( Primitive-Part) ; 
REAL; 

REAL; 


SURROGATE; 

PART-OF(Object); 

Classification; 
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We  discuss  the  following  query  based  on  the  relations  OBJECT  and  PART  described  in  the 
previous  section:  Show  a  material  list  of  the  article  with  number  1200. 


with  implicit  join: 


without  implicit  join: 


SELECT  Material 
FROM  Object. Part 
WHERE  Art#-1200; 


SELECT  Material 
FROM  Object,  Part 
WHERE  Art/(^-1200  AND 

Object .0#“Part.0#: 


The  linear  implicit  join  from  OBJECT  to  PART  is  an  equi-join  between  parent  relation  and 
direct  child  relation.  The  notation  for  implicit  joins  also  generalizes  to  subschemes  which  are 
hierarchical  rather  than  linear.  Precise  definitions  and  illustrative  examples  are  given  in 
[Meier  and  Lorie  1983]. 

An  example  for  using  the  built-in  function  KEY  is  based  on  the  primarily  defined 
KEY-INDEX  combining  the  SURROGATE  attribute  and  the  US^  KEY  Art^  m 
relation  OBJECT.  We  consider  the  following  query:  Give  ail  article  numbers  of  objects 
which  comprise  metal  parts. 

with  KEY- INDEX:  without  KEY-INDEX: 


SELECT  KEY(0#) 

FROM  Part 

WHERE  Mater  1 a1 =metal : 


SELECT  Art# 

FROM  Object,  Part 
WHERE  Material *metal  AND 
Object. 0#»Part.0#: 


It  should  be  noted  that  even  if  the  key  value  is  null,  the  built-in  function  KEY  can  still  be 
performed  since  only  identifiers  are  needed.  Of  course,  the  proposed  KEY  (and  SURR) 
concept  is  minimal  and  helps  to  avoid  writing  additional  joins  to  retrieve  user  keys.  It  does 
not  help  when  more  than  the  'user  key  is  desired  from  referenced  relations. 

The  implicit  join  and  the  built-in  functions  KEY  and  SURR  can  be  used  advantageously  for 
insertion  and  deletion.  Although  updating  through  a  join  is  difficult  in  general,  the  clean 
semantics  of  PART- OF  and  IS  A  structures  allow  using  implicit  joints  in  update  statements. 

We  have  di3cu::::ed  a  general  notion  for  collecting  tuples  from  different  relations  by 
introducing  a  SURROGATE  concept  which  allows  retrieving  and  manipulating  structured 
data  Since  the  system  knows  both  the  structure  and  the  internal  representation  of  the  data 
from  the  system  catalogs,  it  can  optimize  the  implicit  join  accordingly  and  decide  whether  or 
not  to  use  the  specialized  KEY  INDEX.  Accessing  solids  as  interrelated  data  in  a  CAD 
database  directly  instead  of  scanning  throug  different  relations  improves  performance  in 
design  work. 


4.  VECTORS,  MATRICES,  AND  TENSORS 


We  have  argued  that  in  geometric  modeling  it  might  be  interesting  to  store  the  objects  with 
all  their  necessary  geometric  and  topological  information  into  a  database.  Whatever 
representation  is  chosen  for  solids,  describing  points  or  vectors  in  coordinate  space  should  be 
possible.  Also,  since  transformations  are  common  to  all  geometric  modelers,  mappings 
should  be  supported.  The  relational  model  allows  only  atomic  values  as  attribute  elements, 
therefore  it  should  be  augmented  [Meloni  1985]  to  capture  vectors,  matrices,  and  censors. 

A  structured  type  of  rank  m  is  given  by  m  indices,  a  dimension  vector  n^ . nm,  and  each 

element  of  that  type  has  ni*...»nm  coordinate  values.  An  index  is  a  sequence  of  INTEGER 
values  and  ranges  from  1  up  to  its  corresponding  dimension  n)j.  Each  coordinate  is  an 
INTEGER  or  RE.AL  value  and  may  be  un-'  ;uely  identified  by  a  combination  of  index  values 
out  of  ii . i^. 

Examples  of  structured  types  are  given  in  Fig.  4:  A  structured  type  of  rank  1  with  index 

l,...,n  is  a  vector  of  dimension  n.  a  type  of  rank  2  with  indices  1 . n^  and  l,...,n2  is  a 

ni*n2-m.atrix  and  so  on. 


Fig.  4  Dimension  and  Rank  of  Structured  Types. 


There  arc  two  kin-ds  of  operations  for  a  structured  type  U  of  rank  r^  with  dimension  vector 
ni . and  type  V  of  rank  ry  with  dimension  vector  mi,...,mry  respectively.  The  first  class 

of  operations  leads  to  a  result  of  unchanged  rank  and  dimension:  Addition  U  +  V  and 
subtraction  V-U  where  the  precondition  ru  =  rv  and  ni  =  m|  must  hold  for  all  i  from  1  to 
r,j  =  rv:  finally,  nzuitiplicatxn  s»U  and  division  U/s  where  s  is  a  scalar.  The  second  class 
con.  ;,  m  of  3  single  operation  with  a  result  of  indifferent  rank  and  reduced  dimension;  The 
inner  product  U»V  if  the  precondition  ru  =  rv  and  nry  =  mi  holds.  The  new  rank  is  given  by  the 

formula  ri,.v  =  ru-f-rv-2,  and  the  dimension  vector  of  U»V  Is  built  by  dropping  the  last 
com.poncnt  of  the  dimiension  vector  of  U  and  the  first  component  of  the  dimension  vector  of 
V  and  concatenating  the  rest.  As  an  example,  the  inner  product  of  a  n»m-muri.x  U  and  a 
mi»k-ma*rix  V  resulLs  in  the  new  n»k-matrix  U*V. 

Sine  ■  we  have  introduced  a  new  attribute  type,  implications  for  relational  operators  have  to 
tvj  ::'.;u;ed.  Thhe  traditional  mathematical  set  operations  union,  intemection,  and  difference 
are  .:Lil!  porsitie  if  the  relations  are  of  the  same  degree,  i.e.  having  the  .same  number  of 
attribu'es.  Also,  the  Carte;;i3n  product  for  relations  with  structured  types  m.  ^y  be  defined  in 
the  u.sual  way. 

The  p’^ojection  operator,  however,  may  now  be  applied  not  only  to  attributes  but  also  to  a 
.'■truciured  type  itself  to  define  a  new  type  of  reduced  rank  or  dimension.  Thcrefoic,  a 
projection  operator  becomes  a  vector,  matrix  or  tensor  constnicior. 

For  the  selection  operator  Sf.<Relation),  we  have  to  generalize  the  formula  F  .slightly: 
Coristants  in  a  fonnula  may  now  involve  coordinates  of  a  structured  type. 
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The  additional  operators  join  and  division  could  be  defined  by  first  "unnesting"  (compare 
[Schek  and  Scholl  1984])  the  strurtured  types  and  applying  the  usual  relational  operators. 
However,  performance  would  become  a  problem.  Also,  we  don’t  plan  to  support  a  join 
concept  for  structured  types  based  on  coordinates.  We  rather  restrict  join  and  division  for 
structured  types  by  rank  and  dimension  conditions:  Two  relations  with  structured  types  are 
called  join  compatible  if  corresponding  attributes  show  same  rank  and  dimension,  and  if  their 
values  are  drawn  fVom  the  same  domain.  Of  course,  structured  types  involved  in  a  join  may 
first  be  projected  in  order  to  make  the  relations  join  compatible. 
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5.  HRST  RESULTS  AND  CONCLUSIONS 

After  a  decade  of  research  and  development  activity,  relational  database  systems  are  now 
available  as  products.  The  flexibility,  logical  simplicity,  and  mathematical  rigor  of  such 
dat.abase  systems  demonstrate  a  significant  new  approach  to  data  management,  especially  in 
the  b’.usiness  application  environment.  Today,  relational  database  systems  are  also  attracting 
interest  from  users  outside  the  commercial  areas  for  which  such  systems  were  initially 
designed.  In  particular,  the  need  for  efficient  management  of  engineering  and  design  data 
ha.s  triggered  resear  h  on  both  the  requirements  of  such  systems  and  on  extensions  to 
existing  database  systems. 

Some  experimicnts  have  already  been  made  by  extending  System  R  [Lorie  et  al.  1984]  to 
generate  and  support  surrogates  for  engineering  applications.  System  R  catalogs  have  been 
modified  to  capture  the  structure  of  a  complex  object.  The  structure  information  allows  the 
sym.em  to  analyze  the  implicit  join  operator  and  to  find  all  necessary  links  in  order  to 
materialize  the  query.  Also,  a  special  system  table  is  maintained  for  each  hierarchy  of 
relations  to  im.p’em  :nt  a  fast  intra-objecl  access  path.  This  path  is  used  to  enforce 
parent  child  integrity  constraints  and  provides  better  performance  for  clustered  access  and 
m.anipul alien  of  tuples  which  be’  mg  to  the  same  complex  object 

We  have  im.p'.em.entcd  a  3D  modeler  [Meier  et  al.  1985]  based  on  a  hybrid  data  structure, 
namely  'ustng  a  CSG  approach  for  the  user  interface  but  storing  the  designed  objects  in 
boundary  repre.teniation.  Our  .'^olids  are  restricted  to  plane-faced  objects  however  (see  Fig. 
5)  The  u.'^er  .m.ay  translate,  rotate,  or  scale  objects  or  may  choose  a  Boolean  operation  for 
union,  mtemcction,  or  difference.  Also,  a  hidden  line  algorithm  is  available  to  better 
visualize  the  objcci.s.  This  modeler  POLY  has  been  combined  with  a  relational  DBMS 
dev'./icped  at  our  institute  in  order  to  study  interactions  between  geometric  m.odeling  and 
da  tab. as  es. 


Fig.  5;  Using  a  3D  .Modeler  for  Defining  Plane-Shaped  Solids. 


At  present,  we  are  te.'-ting  a  .storage  structure  [Durrer  el  al.  1985]  for  complex  objects  and 
vemicn.^  rased  cn  the  surrogate  model.  A  multiple-tuple- at-the-time  interface  for  PART  OF 
and  IS- A  structures  !.s  the  most  important  distinction  to  conventional  DB.MS.  This 
prac  ;d'ura]  interface  allows  us  to  fetch,  copy  or  delete  a  complex  object  by  a  single  database 
call  ;n  order  to  ran  geometric  and  graphical  applications  more  efficiently. 


Acyr.oy.:,^d'^ru:n:3 .  I  am  grateful  to  Konrad  Durrer,  Erwin  Petry,  Fmwin  Reiner.  Walter 
Reim.  r,  and  Andreas  Walchlin  for  implementing  a  storage  structure  based  on  the  surrogate 
mode!.  Helpful  comments  on  an  earlier  version  of  this  paper  by  Klaus  Hinrichs  and  Hans 
Hintcrberger  are  also  acknowledged. 
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Fignres 

Fig.  1:  Using  Database  Techniques  for  Solid  Modeling. 

Fig.  2:  Conceptual  Scheme  of  the  CSG-Approach. 

Fig.  3;  Conceptual  Scheme  of  the  BR-Approach. 

Fig.  4  Dimension  and  Rank  of  Structured  Types. 

Fig.  5;  Using  a  3D  Modeler  for  Defining  Plane-Shaped  Solids. 
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Fig.  1:  Using  Database  Techniques  for  Solid  Modeling. 
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\  oivlc-r,  !  e  .  X  <x  for  all  i  F-urther  \xc  a'.sunu-  without  loss  of  itciicrality  that  x  <x  !or 
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In  ih's  section  we  study  scxeral  special  cases  nf  the  visibiliiv  problem.  linear  time 
aleorithms  are  proposed  even  when  the  s  s  are  not  sorted  in  the  a  -  and  h  -  (jrder  1  lie  first 
s.ise  IS  the  rnntre!  where  h  =  f).  a  ><).  i  =  1,2.  .n.  The  second  ease  is  the  (C-i.h/e  lomh 

.  ,;if,  where  ^b  o.y  jp.j)  (,,r  each  i.  .i,  =»  A  or  h  =  O  We  assume  the  s 's  are  in  sorted 


Lemma  5  Determining  whether  two  arbitrary  segments  s,.  s^  from  a  set  ol  input  seemeiits 
l.irm  a  visibiliis  pair  requires  ‘.Mn  log  n)  steps  il  the  s^  s  are  mit  si.rteJ  in  the  a,  -  aiul  b 
order 

/Vo  o  Tirst  we  prose  that,  if  an  algorithms  .V  <'f  orsler  <  ('flu  K'g  n!  existed  for  ’.I'.is 
prob'em,  then  the  problem  i  f  determining  whether 

U 

s  an  intersa:  (  u  spins  mto  seseral  disjoint  intervals)  could  also  be  solsed  in  the  same  lime 
In  fact,  gisen  a  set  of  n  inters.tls  ja  .bj.  assign  to  each  ol  them  a  random  coordinate  salue  \  . 

rid  interpret  them  as  sertical  segments  s.  Then  compute  the  maximum  value  of  ail  a,, 

■Ve  min.murn  .alee  \ of  jil  b.  and  the  mr.imiim  and  maximum  values  x,,,,,,,  x,,,^  ol  ail  \ 
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-  ,  .  thus  sob.ir.e  me  union  pro'  iem  for  ia.bi  Howeser.  this  latter  problem  requires 
'..’'n  .  'g  nl  steps,  as  shown  in  Append. x  1  V'v  c  eon-Jasle  that  .V  eanntu  exist 

|m 
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I  Scan  the  sequence  of  all  segments  in  order  of  increasing  values  of  h,,  and  extract  the 
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a  .s  ■:  I  o  M  s  Ml  a:  .ni  i  ■ 
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N('tc  thal  j(  jn\  time  durinj:  the  al^torithm.  all  the  seuments  to  the  left  of  s  are  at  least  as  tall 


Suppose  I  <u<  <1  for  st'me  m>l  are  deleted  durir.c  s'tie  execution  of  the  interior  ol 


I'te  outer  while  loop  Ihen. 


1  a  ,  >  a .  =  a ,  =  a 


2  )  r'l  1  )  =  i.  ,  for  all  k.  1  <  k  <  m  —  1 
3)  ii^)  =  I  for  all  k,  2<k<in 

Lemma  6  ■kssume  i<|  Then  i  and  |  -ee  each  other  if  and  >  niv  P  i  and  |  hee ome  aJiaeent  in 


Pru.if  (Suffieiencs 


'suppose  I  sees  i  Then  a,<niinta.,a  )  for  all  k.  ’.<k<)^  Suppose  i  and  j  do  not  heeoine 
jPM.'iPt  in  T  hen  there  must  he  s(>me  segment  between  ih.en.  s^ .  :<k<i  whieh  is  nnt 

Jeie'ed  before  either  i  or  j 

If  I  and  I  are  deh-ted  at  the  same  time,  then  by  I  )  a.  =  a,.  and  by  2)  and  2))  k  must  aiso  he 
leleied  However  then  i.  =  a,  =  a,,  vvh'ch  is  not  possible. 

If  I  and  J  are  not  deleted  at  the  same  time,  then  let  k  he  ‘ne  seeinent  aojaeetu  t  >  i  .i: 

1'.  he'. er  one  is  deleted  first,  at  the  time  of  this  ileletion.  I  lien  hy  1  )  >  th.e  heiefil  '  !  tltis 
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sue  h  th.il  a ,  >  min  (a  . a  ) 


('use  I.  a,  <  a| 


Let  k  be  the  closest  seitment  to  i,  such  that  \<k<j  and  >  a  .  Such  a  k  must  exist  il  , 
does  nnt  see  j.  In  this  case  all  seements  between  i  and  k  are  shorter  than  i  and  k.  so  k  will  nut 
be  deleted  until  il  becomes  adjacent  to  i.  But  then  since  >  a,,  k  will  be  deleted  at  the  satne 
lime  or  after  i.  and  i  will  never  become  adjacent  icj  j 


Lci  PC  cluscsi  !•^  j  iP.jt  i<k<j  ,ina  .1^  1*  .1  iPc:i  rn  the  s.iiiic  tc.ispnini:.  k  can.poi 

he  dclcicO  hclpic  ;  -.r  i  .imi  )  lanmit  become  Jcljaccii!  lo  i  in 


I  i:c  C!  I'r,  u  1  lies .  ..|  ’:;e  ,i.c<'r  i;  hm  I-  i.ioes  ^iirc:'.'’.  fr"r.  '.nc  '.cmrr.j,  since  ee.ch  secmeti!  s 
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c'.e;.  i  I  e'  !  =  n;;n  h  , ;  ^  :  .,nj  l_  =  max  I  a  .  1  <  i  <  n  h  (See  Fit:  (■)  c'lnsi'Jer  'he 


-•ni!  .ri  t:  e  ';c;n.,i  se.tmeiit'  eio"  n  a  arels  to  fea^h  1. 
en  f  .'i  r.cn.u  secmen’.s  upAarCs  t.i  reach  L. 
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)  file  s'lnhineil  resuit  svi!!  include  ail  '.isdniils  p.u'i 
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m  ;  o't.rse.:  ..i,'.  tiler  scc.'iieiil s  rHtA.-en  and  s  I  hen  either  m  segments  ,ir 


■  '  he  T',  ;i'  ate 
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The  input  secin-n:-.  s.  ,  =  (v  ,j,h.J  for  .i!l  i.  are  sue  h  that  Oe.,  ^ ^  and  hu' 

eaeh  a  =  A  or  h  =  f)  I  he  s 's  are  in  sorted  x.-ureler.  (Se"  Tip,  7) 

Ihe  se  .Titcivs  .o-  .issuiiud  to  be  separated  into  tsso  lists,  T  and  B,  svhieli  are  doubly 
linked  usine  n'li's  ..nd  (  ;  in  increasing’  x. -order  d,,  is  the  initial  element  of  T,  B,,  is  I.he 
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initial  clctnont  of  B  T  consists  of  those  sepmcnts  with  a,  »  A  and  B  consists  of  those 
'cernents  with  b  =  0  Note  that  a  segment  with  a,  »  A  and  b,  =  0  can  be  in  eithei  Usi.  but 
Will  be  assuitteJ  to  he  in  einly  one, 

!'he  jicoriihr;  for  this  case  is  disided  into  two  phases  Durinc  tlic  lost  phase  tin.- 
'■lb:'. ties  ..itii'i'j  seciT.ents  of  I'  w.li  be  obtaincel.  Tins  is  acciiinplished  as  in  the  rooted  case, 
o.  s.arii-.'ic  !r  '-n  lei:  to  rojht  the  'ist  of  the  F  segments  and  deleting  li'.e  !e(t  'oea!  mjs.riun;. 
■In.s  ;'r. 'ce- s; r. c  Me  T  •  seaments  as  an  upside  do'.cn  rocited  case  Heiweser.  reporting  sisib.h- 
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'  it.le  r  I  s  j  nr!  s  i<  I 

.T.'.';.  aivj  the  :a.  est  B-sea:  ient.  s^  with  has  a^^  <  A 
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•  >!' 

•  a  1 

111:'.. 

iri.TC  p'^ljsc  V'^'c  a '■ 

■  \  c  1 

I:  \  n  ■ 

i  1  > 

^  J 

r  :  .•  '  ; 

r.  iho 

iis[ 

1. 

.  ii  I  hi- 

ric^i  o!  V  I I 

If  \ 

I  Ri'i 

»!H  tile; 

■  Jar' 

Pc  .V 

- 

1  ■  ' 

.f  \ 

[  Kl 

•  i 

I  hts 

■s  f'ccaus-j  's  ris(!) 

.•.ne: 

VC' 

. )  .ire  a 

;  ',v  J s  ■ 

c  c  .’P . 

'  c 

r 

■  'J  ; 

Car. 

t  nu'.t 

I'c  taken  hov.cv.' 

A  . 

sc c;v,c r>  1  s 

’  ocCca 

'cna 

S  . 

.r-. 

•,  ■;-e 

iHJor 

■I"  a 

,,i 

‘  :^c  b 

s  of  the  T-sccmc: 

.,,  . 

,hal  ■ 

s,  an.; 

■',  tirt 

aqu. 

.Mr 

‘  / 

s  '  .  . 

the 

:C 

■f;  ■'! 

i.c  .  '•!  tlic  T---cani 

e  fi  t  s 

'  ^0  se.i 

:;nc.l  it 

1  ipcr 

■ .  0 

",  th 

0  '< , 

rctoD  p 

irccedc'  t.bc  Iclt  or 

0  i ! 

the,- 

P  '  aP. 

ja  ,, 

T  i,  . 

•  . 

'  .  .s 

0  •  ‘ » r 

with  c.: ua!  a  s 

.  •>  Pc'!  He' 

.1  • 

■ir  .  .1  [ 

hcAo  a!ro,is 

V  ho”.'"'. 

•  IP-  , .  • 

- 

ar-'.  -a 

rr.  r-sect 

rt-nts  t-or  w  hieh  V  I  1  1 1 ; 

=  \ 

!R,i'  =  o 

n.oi  .l',fin,' 

a'  ..  i 

1 

t  s  a  • 

■■■■! 

t  e  same  lenatn  In 

U  ' 

case  '.'■e  ..n 

,ih  Jsc  r.i 

■■.  i-tc!  c 

taf 

'■  C  p''OP'- 

'..'Cssuig  to  insert  i 

,  a 

nc 

ct.n.  ; 

'non  there  okisin  t 

IP 

■se_,Tietit 

that 

f 

r  'orir. 0 

J':>  I-s*‘ 

jm-.-nt  on  its  -"jht ,  ■  .ini 

,.-ly 

■  .t  ^  f ' '  r 

s.  n-r  1  .n 

.  1"  ;f;o 

.J  r#. 

.tnJ  1 -s'.'.tiiients  ,‘r,'.ni 

ref: 

to  r:o!'. :  * 

-Id  ■  i-.'pa 

\  ; 

'  *  '  t.  .11^ 

Vs  '1 

1  sc  _”".c:r 

:  N  'le  that  if  V  I  1C  i  i 

;s  t 

lelinc-.l  th',: 

I  Al’i  :  . 

1  i'  .  ■■ 

oh  .s  a 

1  - 

cctticru  an,J  if  k  IRii)  .s  unci, .'lined  ’ 

tl'.en  AP'ii 

■.>.|:i  e,. 

!:-o,  rinc 

)  ■  ; 

SlJ 

.h  M..r 

■1  i  b  ,  lie.  a 'asc  •  ,f  oar 

.sumptions 

a  •’  » s  0  .{> 

■  ,  ■  : 

i-  '  < 

•  r  •• 

.■■inc  I -sc 

■an', cut  k  and  H  s,,'a,'r'.c!; 

'  ) 

'ben  K  f'te 

1  C  !  ■ 

1  •  '•  r  1 

■  ■•, 

At 

1  )  rllilst 

he  present  n  tf'.c  list  1 

.  jn,.';  hc-'tsC'  ' 

.■.."7U':P  '  h 

Pv  CIt 

(  > 

1  Al'ti) 

■sc 

se 

^'tienl  1, 

tS 

the  first 

I  -  seitrncnt  fi,r  evheeh  \' 

I  R 

i  1 ,  1  IS  not  i 

let  Klee!  .Hi 

I'-'.citmcn' 

J  with 

a 

>  b  IS 

■neounterccl  in  the  left 

-  rr' 

-rieht  scan 

w  e  h  .1  v  e 

ail 

•  I  'VhOM 

I  k  ,  .ir, 

e  .nulef ined  1  rorn  i he 

del 

h  ”,  1 1 1  o  n  ' '  ^  ' 

t.  IR  s  -ee 

thit  h  <b  <  <b  So  soiimcnt  j  is  seen,  we  simply  sei  AP(i)  to  j  for  all  i  in 

ii,  !,  I,  I  such  till!  h  This  is  ilone  hv  scaiininu  tile  list  h  .b  .  ,b  in  that  order  until 

I  ■  -  V  ^  i  -  'k 

b  >a,.  1  entss.  !s  s.itislied  or  the  list  is  exhausled.  Iluis.  .it  any  nio'neiu  we  base  a  list  L  L 

■I  1  sean'enl'.  wtiokc  h  IR's  are  undctii'.cd  and  wht'se  b  s  are  in  strict!',  iticreasine  k'rder 

'A”e'j  a  H-secaier.;  ;  is  seen,  sse  simple  eonip.ire  its  a-saUie  esiih  the  b-sah.ic  of  the  firs; 

:n  iCe  is:  I  1  '  ii  s  creati  r.  iheii  Al’ii)  set  t"  j  .ind  ;  :s  deleted.  !r. I  1  ,  l.,e 

['I  '  es,  ref'c.;;'  I  >’ i;  e ;  ks .  sC  .s  e  k!:s..'.!.''ii  ]  .i.ad  eoi.-iniie  scal.i'.itiC  Ihai.i'e  the  sc.i  n .  e.  _ 

.r\  T-skCin.  r.t  kv  ,  Oe'':..'.!  V  FR  .alue  -0  .cc.ored  lexeept  ttia:  t.hc  salue  is  ,;ssie:i'ed  t.i  itie 


or,  esn.  i-.i,,-  J  \  1’  !  'e  M  1 


<  '■  e  I,-;.  '"te.'s  ..re  l.'Ur.t,  the  .iiserlion  'h  .1  i'-N..'ertte''.’.  duriTC  the  kctu-a' 


\  -,'s  :  1  e  ■  .k. 


■;  •  eh  .IS  :  oilow  ' 


A  .  •  It- 


•:  :i  1  I  ■  ■ 


UelAeer  I-  .  'O 


Tie:'.  ‘  ;ts  neiahbi  rs  ire  U  'ceme'  t.,  1'  e  'l  lep.at  een 


•i-e:te  '  ci.- 1  .s e-- 11  t-.so  R.-seenien''.  then  tlie  '.  s.h; or. 


■elv'-o'  a  I*  -  se  ch'.e :t.,n  .'  ■'■.'I  :■  p.rt 


te.eted.  it  .'.1  ici-.h  .rs  ire  not  I'.'th  I  -  ■  r  r-  r; 


poir  .J 


.  .'..e.  ■  •  .:-..t  oe..iii  e  h^ 'e  in.i'.  i'e  e . s  o|  j,.,;,  lehctii.  I'l  Case  I‘ 

I'e  ’'  I,:;,  (  I  Si  ii',.; n;  1  .m,;  ps  left  iieieliboi me  B-seemetU  ,  ■.  I'  lioiied  oiiis  wl.e’ 

h  f  1  'I  -i..!  ij'i.ne  ;  ,ir,(l  ,n  ase  I's  'lie  .isibihfe  of  F-seement  1  .irid  If  -  se  .''lie  n  t  j  is  reporle’ 

•  ■-  'len  .1 


F’hase  Iwii 


For  this  .ilgon’ihrr.  it  is  assumed  that  ii.i-,,  .  .i,,  is  in  the  ordering  of  B-  and  T-segnienis  tn 
a  and  h  eolleetiseis  in  inereasing  order  such  :l-..ii 

1  >  If  s  and  s  are  both  T -segments  aitli  h  =  h  a. id  s  is  lo  She  tiebt  of  s  .  ;her. 

k  Sti  U  'fTI  k  r* 

N  sf  T; 

'  I  If  s  apv'  ■-  are  B-seements  with  a  =  a  a'ld  s  is  to  the  rig!',!  ol  s 

o  '  k  '  O'  a 

s  ai 

a)  If  s  IS  a  B-seement  and  s  is  a  1 -seemen;  siien  that  a  -■  b  .  then  m<k.  i  .  a 

k  ■  k,  V  k 

I  -  -■erne n!  preeed.es  a  B-sk  gment  in  the  oidering  when  then  a-  .:nd  b-values  are  eijuai 

Friproeess 
k-  I 

NPL  — ;  \!'!  ,s  a  gueue  .'.nd  is  imtiah/'o!  t'.i  contain  a  kecrne.nt  t  wiio'C  b-'.  ahic  is 

\^!  iik;  \i’l  and  j«»=-M’L  are  insertion  and  deletion  riperat;  ms  on  API 

respe.l;-.  . 
i '  w  h  1 1 s  • ' :  I 

k  -  k  -  : 

1)0  whiie  k  <  n  and  lypeli,^*=)  and  S'fR(i,^)/0 
k  —  i 

f-  'oa  A  !'  1 

1 !  I  .  ;  '  .  '  =  T  ihe n  \ PI  k»ti_ 

eoe  ‘  Iv“e'’^^i  =  B  ihenj<«APl,. 

Do  wbil,,  a,  >b, 

NPi;/-\ 

;<=APL 
f  r.dwnile; 

A  PI  I  i 
fl 

If  k  an  then  |♦AP1,. 

Do  while  b  #  A  1 
AP(j;-<) 

l♦APL 


f  ndvshile 
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fl 

I  ndu  hilc 

IS  tils'  position  th.jt  IS  imnisiliatclv  to  the  neht  of  the  T-seenient  j  when  j  is 
iiserisil  n  the  list  I.  Keli'w  If  AP(j)  =  ()  then  j  is  to  be  insertsJ  into  the  nehimost  positisin  of 
1  1 

In  t..!:..'s  :l;e  I  ,i.',t.ii:i  the  H-ssenieiUs  in  increasing  x, -order  and  maintain  it  as  ,i 
Oc'uh!'-  -link'Ssl  i;si 
k  -  ; 

O"  A!i:;e  m  ' 

1!  1  >  r.  I  )  =  1  then  \.:iae—  else  Value  —a  fl 
Ho  Atiiie  T'.peti  i  T  rind  h  =  Value 

i. 

i;  .t:"  then  !:iser'.  r  to  the  left  of  APir^l 

si  s  iri-ert  i^  at  the  end  o;  1. 


11  I pc!  f.r ’.. ) )  =  ft  then  Report  (i^.nti^))  fl 
I;  f  '.pet  (i^>)=[i  and  \  Tl  then  Report  (Mi.J.ii^)  11 

If  T;.  pe' m  I.  ’  I --- H  anti  f  t  pet  ■  1 1^^  I  )  =  R 
r.  K.,p'  t!  I  Ip),  ni  i.  )  l  fl 
^  k  1 
f-  nC.'.t  nde 

I/''  A':, A'  i  .;ei._^i  =  l’  .intl  .i  =  \'jlue 

1:  i  ;e!  )  I -•  H  th-sn  Report  (i^.n(i^)i  II 

it  1  .  pel  ;  1.  I  ;=  ii  Report  t  ■  (\  to,  1  !l 

I'  r-.  pa':,'.  I  *!■•  ;'i!  intl  ,i  n'lr.  (i  ,a  ) 

'e  >  '  '  'k ' 

tile  n  e'.n  a  1  ;  ;  r. ( i_  i  i  f  I 

Deiet.  p  o  1"  I 

k  -  k  -  1 
f;  nd'.v  Pile 
fndw  hile 


Double  Comb  »ilhuut  a  or  b  sortings 


For  j  sc'^tTiL'nl,  s,,  ihc  set  of  sisible  segments.  I  )  |  I:.))  or  Ij.i)  is  a  sisibiliiv  pair  (  r 
o.ripl-,  !ht  set  of  nodes  sshieh  heeonie  .ivlj.ieeni  to  n-nlc  i  '.ii  the  :ist  sonutirpe  Jiirine  the  s^an 
II  ;s  pvissihie  to  determine  these  seiinienis  without  aetujHs  performinit  i!ie  'eaii.  and  tlv.i' 
usini:  the  a  or  b,  sorlitle’s 

I  el  K\l,;  =  1  ;  i  0,1  IS  a  sisibdiiy  pair  and  j>il  KVli)  o  the  set  of  seentents  on  :'’e  -  j!': 

:  s  a:,-,;-,  .ite  ,isd'  _■  ;  >  s  ,  Ciearl'..  li  R\'fii  is  uetertronesi  tor  e.ivh  s  'iien  .dt  ,  p,;,o 

.  t'e  repottesl  Sapp'ose  now  that  s  is  a  IF-seitmcni  and  eon.sider  t!ie  f.  ho.s  ;v;  two  se 
,;;.:.ses  d  seeme^t'.  Itii)  and  T(i) 

I'e;.  =  Is  I  see!-!  that  s^  is  a  H-seitnient  ssith  i<c)  r.d  a  <a  for  ol  1  d  It 
Oofi  n  sj  =  ni.i  and  a  <a  for  ail  l<i<k  and  if  s  is  a  if-se-.-njent  w;t!i  o  .-  .js'ii  .  i.n 

•i  •  '  si  S  .  ,  .  1  . 

.  .  1  <  ,  <  k  t  ht  r.  ..  e  a 

ioi'  IS  a  way,  I'e'me  i  uidess  s,  happens  to  he  the  last  seamen:,  in  tl'.e  doj's'.e  e'  nt's  H,;, 
■■  ■  's  ,d  ext^eth.  iho-e  se^rnents  s,  w.-«uld  sec  if  there  were  no  f- segment ,  to  hloek  s  ■ 
•  lienee  ft'i)  ,an  be  determit! ed  for  all  Fi-seetnents  Fiy  jppleine  tne  •  ■e.t-.et-sase  .deor- 

:..f.  ;  tie  ist  of  H-',;;.' merits  while  leno.-i  u:  tfie  T-seenicnis 


’'it  -  is  [  _  A, ere  eaeh  IS  a  I-  eemeni  ,ind  s.,  is  t.he  f.rst  f-seeinen;  .siir,  ,.  p,  .,rs, 
p<p  .  .n^yb  amf  if  s,  is  a  T-senment  sucti  that  p<pe  p,..  '  -a  nu  thm 


'I,  o.j'.  he  err, t)i  .  ihcre  are  no  T-seameiils  to  the  neht  if  .  w'’,.i  e  ;■ 

'  ■■•r!."  •,e.,t  s  s  fti)  consists  of  e.xaetiv  t.'iose  T-seeriient  ,  A.h;,,'' 

'•1  .  i  ■  eht ,  if  there  were  no  other  B-se»;menls  to  block  s  .  e  ision 


i'j  'll,  kk'tii  e.in  be  determined  he  soannmp  hoip  fp  .  in,:  I  ,,  iu 
,  „ ''.p.-.r >  lojir  lfr,,-ri  '.eft  to  neht).  and  haltini’  svhen  eutrer  of  tpe  two  seipie  r.  ee  s  inl.rlere 
wuh  s  s  erion  o!  ttie  other  Idle  definition  below  piees  the  first  pair  of  segments,  one  from 
'■a.ii  .eaiicoee  that  f  lo.ks  s  's  eision  on  its  riithl. 

Ill  f  ne  '  e  h  I  s jrh  t h at  !  s'  ^  <  k.  I  <  r,  <  m  Hid  a  a  h  and  il  p,  <" q  then  b  .01 

-  -  a,~  rs,  * '<  a,  . 

elsi  I  ;■  cj  ^p  .  ihen)  h 


19 


I  ||  aiul  ,  ihcn  Ift  h  =  (»  If  =  0  and  P^<q^  ihcn  let  c  =  <).] 

i;  '-■ji.h  a  p.'.;r  then 

R  \  I : )  =  '  _  tj  i  ,  • 

:  ,1.  ■  ;■  ;,r  .niijs  •  .  ms!  tne;;  i’.  ;  i  .i-ul  li;)  d(  t’ut  interfere  '.sith  eaeh  I'ther  .ind 

!n;.  thi'.  ..-e  dettne  ;  e.’-i  '  =  iN.n'.]  It  is  elear  that  te.h)  ean  he  tcend 

'.eannai-  I',;  anel  H'l)  ^;niu!tane  'esl'.'  anJ  that  the  lutmber  of  seemcnls  examined  is  iV' 

■  ..(  ^  ■  ''o.'f-  t'ai.-s  rep 'fted.  Aith.in  a  e«  ns'.,in!.  ilxaetlv  the  same  concept 

..■■•pees  •  :  ...  it'r..ti'',.e  R\;,  :e,r  I-secnier.t  \\  e  sh.al'.  ei\e  beUv.v  a  brief  descnpii'in  .'f  the 
'.ir  ^''mt'  ititt'c  RVo!  fo;  e  ich  ff-  and  r-secmen:  Nh'te  that  wc  should  no!  crmipute 
ij  j  ,  ,  ;;  '  ,  s,  ff.  .,r  [•sejment,  .'ttier'aise  t.’ie  time  taken  ma'.  he  me.re  th.m 

'  a.  1  K  k  ..  r.t.itn  much,  ‘e-.ser  elements  than  tnose  tont-meit  :n  Hiii 

...  I  I  I  \,  tim.',-  needed  f  T  ■;cirf.putinc’  RV'lii  for  ail  .  is  propi  rf.' >nal  to 

■  ■.  .  n,,,;  ,  •■ ..  rv  -  r,.;:;ed  ;ii  .  me  •otal  number  -d  !>-  anil  r-searr.e:  ts  examined 

■.X,_  ■:  ::  .  ,  ,  .  I-  .m.l  P.-secmenis  in  sorted  x -order  from  richt  eft  Dunne 

.V  .1  .  ,  '  t  ;  I  'S'  'ments  i,...  i.  th.il  '.xe  nace  se.inned  in.  ‘  r.  it  'tder  ixith 

's  e  .ri  ;  1  ,  ;  It-s-.  .m,,  nts  j|.i- . ),p  '-xitli  a^  >a  >  >a,  \khen  a  T-se'cment  t 

;  ;■l;..^.XJ'd  aT'd  I,  Is  deleu'd  from  r  and  inc:!.d.’.J  :n  R\(!i  i( 

;  :;■..  j  it  rn'd  for-xjrd  from  i.  aiul  i.,  is  deleted,  irom  It  .in.d  included  in 

K  X  I  i  .1  V-  ;  ...■■■■■'■  t  1  .ippeiided  'o  ifie  list  f  Wlic:’  a  ft -secment  s  is  seen, 

...  le.en.t  is  .ippenile.l  to  the  list  It  (1  or  m  .  x  iiiit  le  sc-.'  F  ten 


t  (ii-ni  ral  i  ase 

In  nils  '  'etiori  -A,  .on.iiJ'.'r  tne  ceneral  problem  and  propose  ,i  method  t..  .Ir.ide  it  into 
siihpr.  'fiien-'  a,  Pi.ifi  ..ir-;  .ioMOie  eomtss  1  lien  sop.c  eaeh  slo-i.-'le  eomb  e  ase  sep-itiilely  and 
eombine  the  r-sults  i.'aelher  fhe  complete  jlitorithm  takes  time  (J(S)  scherc  S  e.  the  total 


iniinf'er  of  se  tnu  n-.s  ..tie  r  : I' ■■ 


o.  ision 


hik- 


If  N  I'  U'  ho  'rn.ill  onou^h.  :hcn  ihc  prohioni  o.in  bo  o.isiK  dui  iod  into  ih'ubio 

O'-fTib  'i.f^;vii"loms,  u'lno  tfio  ;in(J  ICitil's.  The  suhprnbloriis  wii!  ho  luimborod  frnni 

no  l  -  'on;  'I'.ub  'iibpr>iblorT'.  must  be  provided  with  its  1  anil  H  lists  ol  sooments  m  soitod 
'■  1  '.S'!  I  lios"  'r.sis  .;-o  'si.il!  ■  .multaiu'oiisly  in  mio  so.ii;  vd  tho  ontiro  sot  of  s. v: lo.oni s.  m 
'.•Old-.:  \-o,i..  h  -ojo  •.■lit  '  .  Is  onooiiitiorod.  It  is  appoml-o-J  to  tho  M-'.:st  . I  (mm.  the  I -!m! 

ol  1  ti  ll  .ltd  It-'i  ;  .i-f.  It .iird'or  i'rop,,-r'.\  botwoon  I  I'iir  .ir,d  i.fvii 

i  dt'o'  tfio  '.^dft  ot  uithoui  sorlino  :n.i>  ho  ttsod  to  si'ivo  ilie  siibprobloms. 

'0  'tt  ’0  s.  :',Osos-.ir\  to  provide  soitiiuis  b;.  a.  .Old  b,  C'a.'o  must  f'O  taken  to 

..ills'..  !,-,e  ,1' s .. -i;'; ,  o-.s  I'tt  'iio  s  and  b  vs  fiioli  aro  m.ido  by  ;ho  double  comb  aioorithm. 
lb,  ...It  bv  s.i '.  1- ;  I'.'d  i'..  ..siit'c  .1  bu..k.‘l  ssift  to  order  all  ondpoints  i.d  ooual  v -ooordmatos  bs 
: :  o  _  ...  '.rd.i'.i't  it  •:  . .  ■isidonnc  ai!  b 's  of  .i  ■.■■»on  ‘..iti.o  boforo  .my  .i 's  Vho 

■  .  \  o  -.s  ..  r  .•  ropi.Mod  n.  tho  '  i.t.os  iho  hon.'  "tta.  linos  bouitd- 

.  boon  soi',..'d,.  simto  proOsSstne  nay  bo  ro.p.ii’od  to  rontoso 
no  ,1..  .-omplishc-i.l  in  l.noar  tmo  iti  tpo  ni.mbor  of  so'cmonts.  n. 

i  '.'ft  .110  o,p,  ,  .  ;  pr.  ,i-.i.;ris  in  s  iiioh  N  is  cuararfoo  J  to  bo  i.no.ir  in  n, 

1  enini.i  "  !’  ii-  -■  ."f.-.n-.  .iro  '.j  p,...th  'laan  N<-n 

I’’’"'  I;  ■•s'.i  •  .•  si:  .s  o,  ;  1  (,  j^ltiiii-i  for  .il!  s  Si.[  pos'o  I'.c)-  Inon  1,‘iorc  nuisi  bo 

'I  .iisnl  s  •'!.  n  .s  ; ms’.'-tod  m  I  (Mii-l  .iiul  d.o.'o’oo’  ii;  !  (r/i.i  Sinoo  s  vsjs  insi.-rtod  m 

1  'i  I  <  I  ;  ;  ;  tr,'.-::  .  ti'.,.:  f'o  .o.icor  s^.  bjji  ij);,  is  ..  ,1  p.iss.hlo 

n 

b  ni.birl  .  d  d’.o  nt.iomam  :  .ii  ,  ol  lenotlts  is  then  S  is  at  most  ((.'-i-l)ri  sinco  oaoh 
.o.mont  .an  h.,-  eui  u(i  into  at  rnosi  ('-el  piooes 

1:  'f.'.'  imoum  .f  ^"ni.iir.mont  is  abo  boundod  by  t'  ihon  N  is  at  most  (('-flln  I'lto 
''''■'ih  r  oi  oiher  .ocmonis  'viio-o  pr'.'cotions  onto  the  y-axis  ,iro  ..lisjoint  but  .ill  oont.iinotl  in 


(hjt  <if  a  sfj:ment,  s  ,  i>,  defined  tn  he  the  c<'nlainmeni  nunihei  of  s  t'learlv.  s  wiil  m.i  he 
diMded  into  nu're  pieces  than  its  containment  niimher  plus  one 

I'fiere  are  unlorliinateh  .  examples  in  which  N  can  he  t)(n').  as  illusli  aiei.1  in  the  example  'i 
F  iC  ><.  xxhere  there  are  r.  2  lone  segments  and  n  2  sti(>rj  veernents  F  u‘-i  -hort  seemen:  a.l, 
ixe  d'xideO  into  ixxo  -xhi.e  each  htnit  one  iiilo  (n  2+1)  T!ie  total  .s  "  *  n  2  +  i'n  2  =  :f  1 


I  he  dividinc  ilc'oriil'.r.:  presented  at  t''e  hecinninc  of  this  -.ecliori  xchi^ti  lakes  t'nn!  In-.e 
.an  ^e  jseJ  to  Mot  o'lt  \x  ;i,p.  the  vx..ct  '.aluc  of  N  i-.  for  anx  specdic  input  hetote  a.ta.;,'. 
riinninc  the  x' p, '..e  dconthm  If  N  terns  out  to  be  lare.-.  sa>.  larger  than  n  loe  :i.  then  o;-,,. 
m.;x  .xani  t('  use  t'te  fopouirtc  t)ln  n)  time  algorithm  instead 

I'  or  i  i  r  '  n  .i.c  rithm  is  desired,  men  tile  alconthm  ‘sceome  .  cjuitc  simple  s::;,e 
I'l  tore  may  "e  spent  t  i.-tsert  a  sectr.ent  mto  a  list  Tlic  sortinit  <•'  the  a 's  and  h  \  is  -.o; 
.■•sent:.:!  sin.e  a  '.an,  -s',  \  .an  ne  used.  hoxvc'.iT  the  'iiiiple  .ileorilfini  presented  here  Uses 
•'■.est  sor'incs  .o  .xe'l  is  :'te  \  -soittre 


The  Secmsr.is  a  II  he  s.ar.nei!  ir  de..  re.o,  in  it  -soordinate.  a  list  o|  ,;il  seetnents  iP. 
\  orPer  eros'.ifij  a  tt  .o.-op.ta!  s.jp.  one  -eill  he  kept  In  order  to  insert  a  nexx  see.mi'nt  ,rio  :;te 


Tee  a,,i  'T 


eh  r.oi.’i  .>.1.;  have  an  emptx  iiac  .ind  each  edpe  xxi'l  ai'O  *i.ix< 


eirt'i,  t',;e  oh  T.  .’.dl  -'e  sei  i(  .,j:  noi.ies  in  the  siihTee  helovx  ;t,  .ire  empf. 


■step  1  Fteild  a  i  Mhirseil  h.:i.,r\  tre  i.  usme  liie  n  sePrr.ents  o^rdered  hy  x  s,  ,is  ii  'ilcs 


Step  2  t  onsi.ier  ea_''  \  xxh.ieh  .s  either  a  h  or  an 


S I  «■  p  2 . 1  F  o  r  e  a  s  h  a 

Find  .1  neiehhor  in  the  list  and  msert.  by  startini:  at  i  in  the  tree  and  usinit  the  empt'. 
fi.ics  to  find  .1  neii’hhor  f<esct  i  s  empty  fl-nt,  and  the  empty  fl.iits  of  the  ed.es 
between  ■  and  the  root 

Step  2.2  f'or  each  a  «y,  report  xisihihtv  pairs  *iih  ns  neiithbors  in  the  list 

hlep  2.3  For  each  b^«y  report  .isibilitv  pairs  of  its  neiithbors.  if  neither  neichbor  h.is  h  =\ 
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Delete  1  from  the  list  Set  I's  empty  flag  and  if  both  edges  from  i  (below)  arc  empty, 
trase'  Ironi  i  to  the  mot,  setting  empty  flags  until  either  a  non-empty  node  or  an 
adiaeent  non  e!i:pi\  edge  is  encountered. 

Repe.it  drill  .'ill  segment.'-  have  been  inserted  and  deletctl. 


('(inelusiiins 

oi  d  :s  rape:,  a  e  sta  Jieu  a  Msibilits  problem  arisine  from  %'LSl  lavout  compaction 
.i.O"a'  ,  niple'..;'.  ...sdiis  hj\e  been  obtaineel  .-Nssuming  the  input  being  sorted  in  their  \- 
'd  ■  . .  ,11 I'a! es .  a  e  p:  'posed  linear  lime  alevinthnis  fi't  s.irious  special  cases  and  ptcsenleef 
.lie  :  -r  lire  sieneral  r'robiein.  whose  time  complexity  wa.s  linear  in  the  number 

im.ii  '■e'cmen',' 


IImw  .  ;  rein, ■sned.er  tli.-  cencral  ■.  isibiiitv  problem  c.in  be  solved  in  linear 


Appendix  1.  Proof  of  Lemma  5 


Several  prohieiH'  telated  Ic  ihc  atiu'ii  of  k -JiTncnMon.il  inUTv.ilv  hjw  'h-v.'i  raisej  i:i 
vom poMtional  eern'.viix  The  orijiinal  forn'.ul.Uions  can  be  fouml  '.n  [s|.[')|  |n  particuiar,  for 

n 

K  -  i.  K'cc  llic  question  whclhcr  the  moavurc  of  U  [a  .b  .  v.an  be  found  in  !c^s  ihan  C)(:i 

;=  I 

‘  vicps  ['ij  \  11',  L.ilive  jnvAcr  to  Klee''-  vjiiesluui  u.is  pro\ i^icd  bi  I’oe ..iry.oi  ariv!  We;,!.-  i.n 
j!io  .'.ticre  ttiat  the  nieasure  piohlem  riviuire-.  ff  i  r  l"'j  u'  under 

l)oc  ',.in  I.-ve  .or.ip.c.O'onai  tuodel  \vi:h  Linear  C.iniparu-onv  iLll  e  i  In  Li!  (',  c.iL'  . 

■  a  d..,-is:  'n  :!ee  lo  'aneled  '.viih  j  linear  function  of  the  'npiits.  'ahich  :v  compi.o.'J  or, 

r.o.  of  the  nod;  f’!',;'  value  \  of  this  function  indicates  h(. tlier  to  or.K'.eij  l''  tne  s  , 

•  )■  or  :''e  .'i.th:  '■  n  i <(!>,  or  \  is  taken  as  the  .iiiswer  ahen  co-niputcsi  ..t  .i  e.i'  L’ 

U;  n  probiein  ea..h  'e.;!'  is  i.iheied  with  '  yes”  or  no’  in -tead  ot  a  Itineti'”' 

i-'  •‘•o  .irt'en,.  s  s  ,  insider  the  ( 'OM I’.ACTNFSS  Jecisio.n  pr.  bieni.  i'nis 

a  'J  '  a  ,'t’  |  is  an  .nterval.  or  ;t  it  splits  mtw  sever,;'.  Ci',  ■ 

I  iNfss  ...it  ''c  trivially  so.b.ed  in  Ol  n  loc  n)  steps,  by  sor'inc  the  'nte'v 

:  '■  'j'.'t.  .nd  ’'adPin;.:  il’e  interval  unon  by  addintt  one  intcrv.il  at  .■  o.r.; 

vv .  ;  e  I'.v'c  f;..i  '  f  i  ■  problein  a..iujliy  requires  U  (n  loc  n)  steps,  i  )ur  irpv-niv  ".t 

i:  .se  I  :r  ■'''  2  .1  :ne  main  theorem  if  !  |(>) 

(  1.J  m  (  '  AI  f'.'S  I  '  I  N  I  S‘s  requires  '.I  ( ",  loe  n )  steps  in  Dll  (d  rr.ode! 

'  1  I  ;  die  ci:  in's  j  .b  ,if  the  n  input  iiner'...!',  be  tre.ited  as  coord  n.it-.  s  .  jn- 

!  r-  sp:se.  '‘1..S  est.iMishinc  a  one -lo-nne  correspon'de nee  between  irri  '  d  ...'s  ,  |  n 

I',:.  I  ,  .ir'i  ;t.  o'.  •  ,.j,  n  .'i  -pire  Vs  e  will  use  the  Icr.ms  se's  of  .■nte-',  .p  . 

•  :  ii  C;  .."iv  1'  '  poi'its  i’  P  ir.i-ce  the  same  path  v  irt  a  D  I  l,fd  a.ci  r'lhm  'or  i  t;M 

d '■,(  IN!  sS  h.  i.i  ill  ihe  same  leaf,  and  receive  the  same  .inswer  yes  or  n  In 

■I  .  1  I  ear  ns  ,  ..eiis  ■.•.i.oi,  .ler-ed  aionc  y  define  a  convex  recion  L  in  the  -p.m',-  :o  .v  i  'i  1' 
no  p  b'.'i.  nc  :  r  '  .th  all  the  points  of  the  straiciil  IiP't  sccnient  p'p  \d  o' 

1  A  ,  ,  I'o  s.oi'e  .iitswer  from  die  .ilv:..rithm 

(  onsid.er  no.v  .rie  n  intervals  «  |0, 1  |.I ,  «  (  |  ,2 !,  .  .1^  »  1  n  -  l.n],  Tl.ere  ..m  (:i-l) 
.ne'rrui:  Jl-'ons  -  ,r,  ,  ,,  of  such  intervals,  such  that  -(1)  o  [,  for  ail  i  (  I'lit  is,  die 

first  nierv.il  is  jlw.ivs  n  the  first  posiiioni  H.ich  of  the  above  perniu’aiuins  r  corresponds  to 


!!■  we'.er. 
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n/2  LONG  SEGMENTS 


r- 


-A- 


J 


- y - 

n/2  SHORT  SEGMENTS 


Fig.  8  Number  of  final  segments  is  Ofn-) 


t-sfCMENTs  iifle  scudwe  *Jort 
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PBQCESSING  phase  TWC  (WITH  SORTING) 
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Fii:  7  (Pjri  2).  Doubli;  comb  case  ('.viih  soriini:) 
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An  example  to  il'.'.ixir^.lc  Lemma 


TIME  TMaCUjM 

tme  cutes 
^MIEE  LOOP 


LIST 

eCEOSE 

ExrC'JT’GN 


c 

BEEOBE  SEGmENTIS) 

EXECUTION  CElETEO 


1 

2 

3 

4 

« 

> 

6 

7 


8 

7 

6 


jS) 


5 

3 

z 

I 

STOP - 0 


^2 


*6 


7 


54 


^8 

»3 


VISIBILITT  lists 

VIO)  •  1,2,3 
V(  I  )  *0,2 
V12)  •  1,0,3 
V(3)  *2, 4, 8.0, 9 
vi4)  ‘S.s.e 


V(S)  •  6,7,4 
V16)  •  3.7 
Vt7)  *6.8,5 
VIS)  •  7,4, 3. 9 
Vt9)  *9,3 


VISIBILITIES 

REPORTED 

10. 1). (1,2) 

!C,2),  12,3) 

(3,6i,(6,7) 

(7,8),(5,7;,(4 

13.4), (4,8) 

13,81,(8,31 

(0,31,13,9) 


f  I.-  X  I'xjiTinii-  "I  l.'if  rudlml  ease 
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r-» 

u 

S- 

i 

.in  orilcred  set  of  inier\als,  hence  to  a  point  P,.  anti  all  these  points  receive  answer  "yes  "  in 

DTLC,  because  U  I  is  clearly  an  interval.  We  prove  now  that  no  pair  P,.P,  of  such  pt>ints  can 

trace  the  same  path  in  DTt.C  In  fact,  if  this  were  the  ease,  each  point  P  on  P^P^  would 
receive  answer  ’ye.s'  It  is  straightforward  to  note  that  any  such  point  P  corresponds  to  an 
ordered  set  of  intervals  of  unit  I  igth  However,  if  r>  1  is  the  smallest  inde.x  value  for  which 
rtki  =  Ij.r  (k)  =  1  .r<s.  for  some  k,  as  soon  as  a  point  P  starts  slit'ing  from  P^  to  P  .  the 
union  of  the  inteiv.ils  of  P  breaks  up  into  two  intervals,  around  the  ci  ordinate  value  r-1. 
Hence,  ilierc  is  at  least  one  point  P  on  P.P^  whose  answer  must  he  iio  against  the  hypothe¬ 
sis.  that  IS,  P,  and  P^  must  correspond  to  different  leaves 

's'*  e  sonclude  th.at  LJT  L.C  must  have  at  least  (n-1)  distinct  leaves,  that  is,  us  depth  is  1.'  (n 

log  n  1 


lac  jhuve  proi  f  makes  use  ol  segments  of  equal  length.  Therefore,  it  ,iIso  proves  that 
)M  r  \(.'TN  TSS  rem.iins  ot  '.f  in  log  n)  in  the  restricted  case  a,  -  b,  =  c,  l<i<n  (i.c,,  all 


Stereophenomenology 

stereo  -  solid 

phenomenology 

-  theory  of  representation 


Representing,  manipulating  and 
reasoning  about  physical  objects 
electronically. 


What  does  science  include? 


representation  of  objects 
surfaces  and  solids 
functional  dependency 
hierarchical  view 
abstract  models  -  features,  etc 
generic  objects 
internal  structure 
flexible  and  nonrigid  objects 

algorithms 
display 
intersection 
motion  planning 

user  interfaces 
editing 

interactive  graphics 

attribute  grammars  -  simplifying  local  changes 

reasoning  about  objects 
grip  positions 
external  forces 
shape 

design  for  function 

manipulation 

gripping  strategies 
object  motion 


pipefradius,  thickness,  length)  :=  (cyll  -  cylS)  H  III  H  112  where 
begin 

cylt  :=  ycylinder(radius-hthickness); 
cyl2  :=  ycylinder (radius); 

III  :={y>  0}; 

=  {y  <  length}; 

=  (cyll  -  cyl2)  H  {y  =  length); 
bottom  :=  (cyll  -  cyl2)  Cl  {y  =  0}; 
outside  :=  cyll  .surface  D  Hi  D  H2; 
top.in_cdge  top  H  cyl2; 
bottom. in_cdgc.  :==  bottom  f)  cyl2 
end; 


112 

top 


jla7igf  d_j)i]>(  ^radius,  thickness,  length: 


or;  oof h_r/ larh (nipple. top. ui_(  dye,  sjlar.yr  .u  :o.:  jm.in^cayc ,  f 


negin 


:=  !;  -frii rhncss+hole_c:e''''i;  tv-; 


jlo'iye  jiipc (radius,  ft.  f  j: 


sjlanye  srnoothfflange. outside,  Jlarigr.hcKo  oi,  ij  !(>:: 
nipple  ;==  pipe(radius,  thickness, 
bottom  :=  ri2j>plc. bottom; 


Top  :=  yP o n n(  .top 


end; 


ovoid(rl,  r2,  rS)  ;=  smooth(cyll  H  cyl2  f]  eylS,  rS)  where 
begin 

cyll  :=  y7novc(-(4  +  9l  16),  zcylinder(r2)); 
cyl2  :=  yrrLove:(/'+Oj  16,  zcylindcr(r2)); 
cyl3  :=  zcylinder(rl); 


x-slice{width)  :=  HinH2  where 
begin 

Hl:  =  {x  >  0}; 

H2:  =  {x  <  width}; 
left:  =  {x  =  0}; 
right:  =  {x  =  width} 
end; 

y-slice(length)  :=  HinH2  where 
begin 

Hl:  =  {y>0}; 

H2:  =  {y<  length}; 
front:  =  H2; 
back:  =  Hi 
end; 

z-slice(height)  :=  HinH2  where 
begin 

Hl:  =  {z  >  0}; 

H2:  =  {z  <  height}; 
top;  =  H1; 
bottom:  =  H2 
end; 

cuboidf  length,  height,  width)  :  = 

x-slice( width)  fl  y-slice(length)  Pi  z-slice(height)  where 
begin 

front. right-edge  :=  front  n  top; 

front. right. top- vertex  :=  front  fl  right  D  top 


vertexl:  =  (x-,  y-,  z-coordinate); 

vertexS;  =  (x-,  y-,  z-coordinate); 
edgel:  =  line(vertexl,  vertex2); 

edgel2:  =  line(vertex7,  vertex8); 

front-face;  =  patch( edge  1,  edge2,  edgeS,  edge4); 

right-face:  =  patch(edge4,  edged,  edgelO,  edgell 


Comparison  of  techniques 


abstract  models  versus  solid  models 

solid  models  versus  surface  representations 

polygonal  patches  versus  bicubic  patches 

versus  algebraic  surfaces 

numerical  techniques  versus  symbolic 

More  general  objects 

nonrigid 

plastic  flowing  into  mold 
sail 

shape  determined  by  external  forces 
coil  spring 

generic  or  parameterized  objects 

Design  for  functionality 

can  objects  be  represented  by  function 
rather  than  shape  and  dimension 


Theorem  (simplest  form).  If  there  exists  a  motion 
of  two  objects  from  an  initial  position  where  they 
are  in  contact  to  a  final  position  where  they  are  in 
contact  then  there  exists  a  motion  whereby  the 
objects  remain  in  contact  at  all  times. 

Motion  is  continuous  but  point  of  contact  is  not. 

Motivation 

1.  Trying  to  place  problems  in  PSPACEMinkage 
motion,  block  motion 

Simple  multiple  object  motion  planning 


3.  Compliant  motion 


HOW  DO  WE  REPRESENT  GENERIC  OBJECTS?  POSSIBLY 
BY  REPRESENTING  ABSTRACT  OBJECTS  WHOSE  POSITION.  SIZE 
AND  SHAPE  ARE  INSTANTIATED  ONLY  WHEN  SPECIFIC  INSTANCE 
NEEDED. 

ALLOWS  MOTION  PLANNING  AND  OTHER  ALGORITHMIC 
TASKS  TO  BE  CARRIED  OUT  FOR  CLASSES  OF  OBJECTS  RATHER 
THAN  INDIVIDUAL  OBJECTS. 


AN  OBJECT  IS  A  PARAMETERIZED  MAP  FROM  A 
CANONICAL  REGION  OF  R^  TO  R^ 


MOTION  IS  A  CONTINUOUS  MAPPING  FROM  10, 11  TO 
PARAMETER  SPACE. 

EXAMPLES 

TRANSLATION 

ROTATION 

GROWTH 

CONTINUOUS  DEFORMATION 
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Lemma:  Configuration  space  is  path  connected. 

Ic 

Lemma:  Hq(S)  =  Z  if  k  path  connected  components 
Lemma:  H^(S)  —  <j)  space  contractible  to  a  point 


Mayer-Vetoris  Theorem: 


/i  1 

Hq{A  HB)-^ 


is  an  exact  sequence. 


^3 


exact  sequence  image  of  h-  =  kernel  of  h- ,  . 


conditions  for  theorem  to  hold 

1.  Translation  is  needed  in  order  that  configuration 
space  object  be  path  connected. 

Alternative  Lemma:  There  exists  two  paths, 
one  in  free  space,  one  in  the  configuration 
space  object. 


2.  Space  must  be  contractible  to  a  point. 
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